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HE FRICTIONLESS STATE OF AGGREGATION 


By K. MENDELSSOHN, 
Clarendon Laboratory, Oxford 


MS. received 26 March 1945 


BSTRACT. The analogy of superconductivity and liquid helium II has been analysed 
| a purely empirical basis. Frictionless transport is considered to be the fundamental 
enomenon. It has been concluded that the ‘‘ superconductive’ electrons and the 
superfluid ”’ helium atoms form an aggregate in momentum space of zero thermal energy 
-state). ‘he z-particles occupy a set of lowest quantum states which is separated from 
e thermal states by an energy gap ; they have zero entropy, but an appreciable zero- 
int energy. The observed specific-heat anomalies are due to a lifting of z-particles 
er this gap. ‘The transport phenomena are caused by a diffusion of z-particles under 
eir zero-point momentum. ‘This zero-point diffusion does not result in dissipation of 
omentum. It explains why frictionless transport is independent of the strength of 
ternal forces and explains also the critical transport rates (threshold value and film 
insfer). Bose-Einstein condensation fails to account for these characteristic features 
cause it does not lead to a zero-point energy of the condensed phase. A number of 
periments to test the hypothesis are suggested. 


SIN RO DUCT LON 


“™\ OME years ago J. G. Daunt and the author (1942) drew attention to the 
strong similarity between the phenomena of superconductivity and liquid 
helium II, which becomes apparent when we consider quite generally 

e properties of particle transport in these substances. In both cases we 

counter frictionless transport which obeys similar rules in spite of the fact 

at the particles are atoms in one case and electrons in the other. his, together 
th an equally striking similarity in the entropy functions (Daunt and Mendels- 
hn, 1945) has led us to advance the hypothesis that both phenomena are 
used by much the same mechanism. We have described this mechanism as 

e appearance, at finite temperature, of particles which are energetically at 

solute zero and do not exchange momentum with the remainder of the substance. 

view of the peculiar way in which the thermally unexcited particles are separated 
ym those in the higher quantum states—they do not form a separate phase— 

» have described the assembly of these particles as an aggregation. 

It is the aim of this paper to discuss in more general terms * the hypothesis 
the basic identity of superconductivity and liquid helium, and the properties 
the aggregation of particles at the lowest quantum levels. It must be emphasized 
it these considerations are not based on any particular theoretical conception 


* A detailed discussion of the energy state of the electrons in superconductors has recently 
n given elsewhere (Daunt and Mendelssohn, 1945) and a detailed account of the corresponding 
snomena in liquid helium will be published later. 
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] 
but that an attempt is made to deduce the properties of the frictionless state 
empirically on the basis of the existing experimental evidence. | 


§2. FRICTIONLESS TRANSPORT | 


Let us first consider the position as regards the experimental facts.* In 
the case of liquid helium the last years before the war had brought to light a 
considerable number of seemingly unconnected phenomena, the most striking 
of which are : extremely low viscosity, very high heat conduction and, in some 
way related to the latter, a thermo-mechanical effect (sometimes called the 
“fountain effect”) which constitutes a flow of helium in the direction of the 
higher temperature. In addition there are the transport phenomena in a film 
of thickness about 5 x 10-® cm. which covers all solid surfaces in contact with 
liquid helium II. The somewhat unsatisfactory state of the experimental 
material is due to the fact that internal friction and heat transport in liquid 
helium II are evidently very complicated processes. ‘The efforts to determine 
such properties as viscosity and heat conduction have consequently led to complex, 
and in many cases contradictory, results. On the other hand, the film transfer 
which at first sight appeared to complicate matters has proved to be of great help 
since the phenomena in films seem to present a much simplified picture of the 
transport effects in the bulk liquid and have the added advantage of yielding well 
reproducible results. 

With regard to superconductivity, the position is more satisfactory on account 
of the greater wealth of confirmed experimental results. Since 1933, when it 
was discovered that zero resistance is coupled with zero magnetic induction, 
there has been a tendency to describe superconductivity as a purely magnetic 
phenomenon, some authors going so far as to consider the disappearance of 
resistance to be due to complete diamagnetic polarization. Such an explanation 
is, however, hardly tenable in view of the recent determination of the Landé 
factor by Kikoin and Goobar (1940), which shows that the super-currents are 
caused by movement of free electrons and are not merely simulated by electron 
spin. We stress this point since it will appear from our discussion that super- 
conductivity is primarily an electric and not a magnetic phenomenon. The 
connection between electric and magnetic effects in superconductors has been 
expressed in rigorous form by the work of F. and H. London (1935). Their 
fundamental equation which connects the magnetic field with the current density 
does not, however, imply, as has sometimes been assumed (Burton et al., 1940, 
p. 320), that the magnetic effects are the primary ones. The London theory is a 
phenomenological one which does not differentiate between cause and effect. 

We have already given a short summary of the corresponding phenomena 
in superconductivity and liquid helium (Daunt and Mendelssohn, 1942), and 
it is not intended to discuss the analogy here in greater detail. Our hypothesis 
implies that the fundamental phenomenon in both cases is that of frictionless 
transport, and we shall now consider the experimental justification for this state- 
ment. In superconductors the evidence is provided by the classical current- 


* Summaries of the phenomena in liquid helium have been given by Darrow (1940) and Jones 


(1940). A description of superconductivity has been given by Shoenberg (1938 d 
of recent work by Jackson (1940). a ea 
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otential measurement (figure 1a) in which part of an electric circuit is super- 
onductive. An ammeter in the circuit indicates a flow of electrons through the 
uperconductor while a voltmeter shows zero potential between any two places 
f the superconductive part of the circuit.* Thus we have a transport of par- 
icles—electrons—through the superconductive metal in which no energy is 
issipated. ‘The corresponding effect in liquid helium II is illustrated by the 
ransfer of helium atoms in a film over the wall of a beaker (Daunt and Mendels- 
ohn, 1939 a) (figure 1b). If the level inside the beaker is higher than outside, 
lium atoms are transferred out of the beaker until the difference in level has 
isappeared. ‘The most remarkable fact about this transfer is that in contra- 
listinction to the action of a syphon the rate of transfer is independent of the 
ength of the path and of the difference in the height of the levels. This goes 


AMMETER | 


VOLT|METER 


YOLONGNOODYIDNS 


a 
‘igure 1. Experiments showing frictionless transport, (a) in a superconductor, and (6) in a 
film of liquid helium II. A sudden change in the level difference does not affect the rate 

of transfer. 


o far that when the level difference is suddenly changed, no change at all is 
ybserved in the rate of transfer. Here again we have a transport of particles— 
n this case helium atoms—without friction, for friction would make the rate of 
ransfer dependent both on the difference in potential and on the length of the 
yath which the atoms have to traverse. [he same phenomenon has been observed 
n flow of the bulk liquid through capillaries (Allen and Misener, 1939 ; Johns, 
Nilhelm and Smith, 1939) although not so clearly defined because it is partly 
nasked by viscous flow of the liquid as a whole. However, in very fine capillaries 
~10~> cm.) the result is practically identical with the observations on films, and 
or wider capillaries it can be expressed in the form 


FoR nAP ee Po) (Sr ohe (1) 


there r is the total rate of transport (volume/time), R the rate due to frictionless 
ow and nAP the amount transported by viscous flow under the pressure difference 
.P. ‘The value of R and its dependence on temperature in the finest capillaries 
ave been found to be of the same order as in films. ‘his means that, as in super- 
onductors, we are faced with the fact of frictionless transport through the rest 
f the substance. ‘There is some evidence that both in superconductors and in 
quid helium the frictionless transport shows a preference for the geometric 


* The most accurate determinations (using persistent currents) show that the resistivity of a 


iperconductor is less than 10~?° ohm-cm, 
25-2 
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surface of the substance—a phenomenon which will be discussed later, and which 
does not affect our general argument. 

The experiments thus show that there are “ superconductive ” electrons and 

SSUpSaiaUs! ” helium atoms which move through their respective “container 
substances” without friction. When using such terminology we must, however, 
keep in mind that we cannot distinguish between individual electrons or individual 
helium atoms, and that these terms can only be used in allotting a certain number 
of particles to a given set of energy states. In the present case the particular set, 
of states is distinguished by the fact that the electrons or helium atoms in this set-— 
we will call them z-particles for short—do not dissipate momentum. 

It must be pointed out straight away that the z-particles do not form a separate 
phase in either case. Apart from the lack of any direct experimental evidence. 
to this effect, the thermodynamic functions, which we shall discuss later, show 
clearly that the appearance of the z-particles is not accompanied by a phase change: 
in coordinate space. It has even been impossible to detect any increase in spatial 
order to account for the observed decrease in entropy which is linked with the) 
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Figure 2. The peculiar position of the “‘ aggregate’ of z-particles ; they exhibit no interaction 
with the rest of the substance, but strong interaction with an external force. 


appearance of the z-particles. There is, however, another way in which th 
assembly of z-particles is distinguished from any other assembly of spatiall 
disordered, freely mobile particles belonging to a particular set of states. This 
is shown by figure 2, which illustrates diagrammatically the relation betwee 
the z-particles, the ‘“‘container substance”’ and an external force (magnetic field,| 
gravitation) acting on* the z-particles. There exists no interaction betwee 
the z-particles and the rest of the substance, but a strong interaction betwee 
them and the external force, which means that the z-particles react to it as al 
separate aggregate. More than this, the z-particles assume under the influence 
of such an external force the aspect of a spatially ordered aggregate whose degree 
of order can become manifest in macroscopic dimensions (e.g. the magnetic fiel 
of a persistent current). It seems therefore convenient to describe the z-particle 
aS forming a particular state of aggregation—since we wish to avoid the term 
phase—which we will call the frictionless state of aggregation (z-state). | 

It must be stressed that this term does not of course refer to the supercon- 
ductive metal or liquid helium as a whole but only to the SSS TED of z- -particles, 
1.e. the “superconductive” electrons and the “superfluid” helium atoms. As 
the temperature of the substance is raised from absolute zero, the z-state weakens 


* As will be seen in §5 the conception of an “ accelerating” force has to be qualified, since 
it seems possible that all frictionless transport is due to zero-point momentum, the acting force 
providing space order rather than acceleration, 
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ind finally disappears. We will show in § 3 that this phenomenon can be 
described as a gradual decrease in the number of z-particles with increasing 
remperature. ‘I'he frictionless state shrinks until, at the transition point or 
he A-point respectively, the last z-particle disappears. In superconductors an 
ndication of the change in population of the z-state is given by the threshold 
Surve which measures the number of electrons that can be transported without 
riction in a given time over a superconductive surface of given dimensions at 
lifferent temperatures. -We have determined a similar quantity for the frictionless 
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Figure 3. Critical transport rates (a) in liquid helium, and (6) in a superconductor (tin). 


ransport of helium atoms in films (Daunt and Mendelssohn, 1939a). It is 
significant that this “‘ rate of transfer’’—the number of helium atoms transported 
without friction in unit time over a helium surface of given dimensions at different 
emperatures—was found to depend on temperature only. Figure 3 gives a 
ypical threshold curve and the rate of transfer in comparison. These functions, 
vhich merely give the amount of z-particles collected in unit time, contain of 
‘ourse an unknown average velocity. If this velocity is independent of the 
emperature, as appears likely (see § 5), the curves can be regarded as indicating 
lirectly the relative number of z-particles as a function of the temperature. 


§3. THERMAL ENERGY OF THE FRICTIONLESS STATE 


The gradual disappearance of the frictionless state of aggregation is accom- 
yanied both in superconductors and in liquid helium II by an excess specific 
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Figure 4. Specific heat (a) of liquid helium, and (0) of the free electrons in tin. 


eat (figure 4). In superconductors it has been possible to separate the specific 
eat which is due to the electrons from that caused by the lattice vibrations 
Daunt, Horseman and Mendelssohn, 1939). It has sometimes been assumed 


| 

| 
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that this excess specific heat, which is caused somehow by those particles parti- 
cipating in the anomalous state, must imply heat exchange between a system of | 
these particles and the rest of the substance. The point has been discussed in | 
the case of superconductors (Burton e¢ al., 1940, p. 298). The mechanism | 
sketched in figure 2 on the other hand wales it very difficult to account for such | | 
an interaction. In fact, our definition of the z-particles—and this definition was | 
derived from experimental results—would not allow direct thermal interaction | 
with the “‘ container substance’ unless we provide for it by additional assumptions. | 
Any such assumption, allowing energy exchange in one case and prohibiting it in| | 
another, would moreover be highly artificial. 

This obvious discrepancy requires an explanation and it raises the question’ 
of the energy state of the z-particles, i.e. of the heat content of the z-state. The | 
specific heats given in figure 4 refer of course to the contributions by the z-particles | 
and by the normal particles, and it is clear that information on the thermal state | 
of the z-particles by themselves must be obtained by experiments of a different | 
character. | 

Such experiments have been carried out both for a superconductor and for. 
liquid helium II. In the superconductor this is the determination of the Thomson ] 
heat (Daunt and Mendelssohn, 1938) ; the experiment and the result have been } 
discussed in detail elsewhere (Daunt and Mendelssohn, 1945). ‘The fact that} 
the ''homson heat of a persistent current is zero means that the superconductive |f 
electrons do not exchange thermal energy with the rest of the substance. ‘This|} 
leads to a peculiar consequence when we consider what must happen in a ring of 
superconductive metal with a persistent current running in it, which is watmed 
up from absolute zero. Since the electrons in the current are not in heat exchange 
with the lattice, they must remain energetically at absolute zero although the 
substance as a whole warms up. 

An idea as to the thermal state of the superfluid helium atoms can be obtained |f 
from any experiment which causes their separation in space from the liquid as al] 
whole. We have seen in the last section that, while there exists no space order] 
generally, an external force has such an ordering influence (see figure 2). In]] 
determining the ‘Thomson heat in a superconductor we have therefore made use|} 
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Figure 5. ‘The mechano-caloric effect in which a separation of z-particles occurs. A flow. 
of liquid helium If through a powder “‘ filter’ is accompanied with a Reetie of the liquid | 
at the entrance and a cooling at the exit. 


of the influence of the magnetic field, and in the case of helium we can use gravi-| 
tation as the acting force. In the experiment shown in figure 1 superfluid) 
helium atoms were separated by surface flow. If the heat content of the superfluid 
helium atoms, like that of the superconductive electrons, is smaller than that. 
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of the normal particles, the liquid drained off must be colder than that inside the 
beaker. We have carried out an experiment of this type (Daunt and Mendels- 
sohn, 1939b), using for practical reasons the flow of liquid helium II through 
tightly packed fine powder instead of flow over an exposed surface. A diagram- 
matical sketch of the method is given in figure 5. As was expected, a heating 
Was observed at the entrance of the liquid into the “filter” and a cooling at the 
exit. H. London (1939) has since given a general thermodynamical interpre- 
tation of this mechano-caloric effect. Under the assumption of complete rever- 
sibility he obtains the following relation’ between the pressure (pf) under which 
the liquid is forced through the power plug and the entropy difference (AS) 
between the helium flowing without friction and the bulk liquid : 


d: 

“o Ce Sn Cree eee (2) 
‘where p is the density of the liquid. Our experiments (which had to be suspended 
in 1939) were of a preliminary character only, but the (dp/dT)-values obtained * 
) yielded a large entropy change indicating that the superfluid helium atoms seemed 
‘to have lost all thermal energy. Fortunately the experiment has since been 
: repeated by Kapitza (1941), using a narrow gap instead of the powder plug. He 
could not only confirm the existence of the mechano-calorific effect and demon- 
/strate with greater accuracy the complete loss of thermal energy suffered by the 
superfluid atoms, but by measuring independently the heat change connected 
| with the effect, he has shown that the process is reversible. 

| It has to be admitted that we cannot yet establish the complete loss of thermal 
energy with the same accuracy for the superfluid helium atoms as in the case of 
»the superconductive electrons.{ Particularly the specific heat at low tempera- 
tures (Pickard and Simon, 1945) has been interpreted as corresponding to a 
“lattice entropy” (cf. H. London, 1939) while frictionless transport has been 
,assumed at the same time. We will return to this question, which is not at all 
clear, in §4. For the present we shall neglect this hypothetical residual entropy 
which, as the experiments show, does not interfere with frictionless transport. 

, ‘The experimental result that in a “‘container substance” which has a finite 
temperature, the z-particles remain energetically at absolute zero but have 
‘momentum, leads to the simple but somewhat surprising conclusion that they 
‘must move without friction. For friction is the dissipation of momentum and 
any particle engaged in a process involving friction must therefore suffer thermal 
excitation. ‘This means that the observation that z-particles remain thermally 
unexcited while being capable of translatory motion is in itself sufficient to 


* The work will be published in detail later. 

+ Allen and Reekie (1939 a) have obtained (dp/dT)-values from the thermo-mechanical (fountain-) 
effect which is the thermodynamic reverse of the mechano-caloric effect. These values are of the 
‘same order as those from the mechano-caloric effect, but again not accurate enough to determine 

the exact entropy loss. 

t Our conception of the properties of the z-state suggests that when a rotating mass of liquid 

helium is cooled below the A-point, a frictionless current of z-particles will be set up which will 
-persist when the momentum of the thermally excited particles has been dissipated. The decay 
of this persistent helium current will provide an equally sensitive measure for the frictional 
resistance of the superfluid atoms as its electric counterpart does for the resistances of the flow of 
superconductive electrons. 
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explain the phenomenon of frictionless transport. It seems permissible to speak | 
of the entropy of the z-state since this state, as will be seen, forms a separate phase — 
in momentum space, and it is clear that this entropy could not retain the value | 
zero—as we know it does—when the state is agitated, unless this agitation does not 
involve friction. . | 
The realization that the entropy of the frictionless state is zero at once removes }f 
the discrepancy with which we were faced when considering the excess specific | 
heat of the z-particles. It also answers the question as to the significance of the | 
entropy change accompanying the weakening of the frictionless state with | 
increasing temperature. ‘There exists no direct thermal interaction between 
the z-particles and the rest of the substance—the specific heat of the z-state 18 | 
zero. The observed specific heat is clearly due to an excitation process in which ]f 
z-particles are lifted into thermally excited states. When this occurs, they cease 
being z-particles and cease to be available for frictionless transport. 
The anomaly in the specific heat connected with this change (figure 4) exhibits 
the familiar pattern associated with ‘cooperative’? phenomena like ferro- | 
magnetism, rochelle-electricity and order-disorder transformations. ‘The latter |}} 
in particular have been made the subject of detailed statistical treatment, the |} 
methods of which have been applied to other phenomena, as for instance molecular 
rotation in solids. However, no general statistical treatment has yet been | 
developed which can be rigorously applied to all these changes without regard | 
to the particular physical process involved. Simon (1930) who first investigated | 
these ‘quantum jumps” in solids has generalized them as excitations of internal | 
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Figure 6. (a) Degree of order, and (b) excess specific heat of different alloys (theoretical). |) 
(1) ‘Type AB with short range order, (2) type AB with long range order, and (3) type AB; |) 
with long range order. Hp is the total energy of the transformation required to change | H 


the alloy from complete order to complete disorder. (Cf. Fowler (1936); Nix and | 1 
Shockley (1938).) i 


degrees of freedom. ‘The large amount of experimental data gathered in thelf} 
last twenty years exhibits all stages of specific-heat anomalies ranging from the} 


uninfluenced Schottky curve over A-points to transitions involving a latent) 
heat. ‘This fact is important for our particular problem, since it shows that no! 
clear-cut distinction is possible with regard to the mechanisms of the various| 
excitation processes, and that to describe an unknown physical process on the}} 


strength of the observed specific-heat function as an order-disorder transformation 
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or a rotational change has little meaning. It usually means only that the par- 
ticular statistical treatment can be applied, but it discloses nothing with regard 
to the nature of the process. In the case of the order-disorder transformations 
it has been possible to connect the theoretically determined states of order with 
specific-heat functions, which approach the observed ones, for a non-cooperative 
(short-range order, AB) type,* a cooperative (long-range order, AB) type, and 
a highly cooperative (AB;) type. Figure 6 illustrates the point that the type of 
the change is not indicative of its nature. The degree of order—the inverse of 
the entropy—is a function, and in certain cases a direct measure, of the population 
of the lower state. 

Considering the apparent lack of space order and the high symmetry of both 
the electron and the helium atom, we are forced to assume that the degrees of 
freedom lost by the assembly of z-particles are those of translatory energy. ‘There 
is nothing unusual in the fact that at finite temperatures some particles in a sub- 
stance are energetically at absolute zero. This is actually true for a proportion 
of the atoms in any solid at low enough temperatures. However, in the case of 
the z-state we deal with freely mobile particles. ‘The nature of the phenomenon 
of frictionless transport is therefore the manifestation of an effect predicted by 
Nernst (1926) when he extended the third law to systems with freely mobile 
particles, i.e. gas degeneracy. 


§4. THE ENERGY SPECTRUM 


While thus frictionless transport is connected with gas degeneracy, this 
does not mean that it is a necessary consequence of it. In fact, the free electrons 
in the Sommerfeld model are in a highly degenerate state without exhibiting 
frictionless transport. ‘The existence of z-particles, capable of translatory motion 
but retaining zero entropy, requires therefore not only degeneracy but degeneracy 
of a particular type. The particular qualification necessary for frictionless 
transport becomes apparent when we consider the term system of the Sommerfeld 
model (figure 7.a). A continuum of energy states is filled at absolute zero up to 
the level (L) which divides a completely filled set of states (A) from a completely 
unoccupied one (B). At finite temperatures an accelerated electron can and will 
dissipate energy, since it can be lifted to any place above (L). If we wish to 
allow for frictionless transport, we must prohibit states bordering (L), so that no 
dissipation of momentum can take place ; i.e. we must introduce a gap in the 
term system (figure 7b). As we have pointed out elsewhere (1945), in super- 
conductors this gap forms the limit of the Fermi distribution at absolute zero. 
It follows from the considerations discussed in the preceding sections that such a 
term system in which the continua (A) and (B) are separated by a gap in the 
energy spectrum must be characteristic not only of a superconductor but of 
any other z-state. (The only other case known is, of course, liquid helium II.) 

The existence of the gap signifies a separation in momentum space between 
the particles in (A) which cannot dissipate momentum and those in (B) which 
are thermally excited and acceleration of which results in friction since they can 


* This function, a Schottky curve, has been observed by Simon (1926) in a number of pure 
metals ; the nature of the change is unknown. Other examples of this type of anomaly are the 
paramagnetic salts and solid ortho-hydrogen. 
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be lifted to any higher state in the continuum (B). We can therefore speak of 
a “condensation” in momentum space, the z-particles forming the “‘ condensed” 
phase. ‘This phase forms the frictionless state of aggregation. Neither in 
superconductors nor in liquid helium is this a condensation to zero energy in the 
model sketched above. This means that neither case corresponds to a Bose- 
Einstein condensation since both have an appreciable zero-point energy. As 
mentioned in §3, it is by no means clear what happens to the “‘ condensed” 
helium atoms in the liquid. The assumption that the rise in the total entropy 
of liquid helium below 1° k., which Pickard and Simon found to be proportional 
to T3, is due to “lattice vibrations”? can be brought in accordance, though with 
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Figure 7. ‘Term systems (a) of the electrons in a normal metal, and (6) as suggested for the 
electrons in a super-conductor or the atoms in liquid helium. A is the zero-point energy. 


some difficulty only, with the phenomenon of frictionless transport, since 1t may 
be argued that the vibrations excited at these temperatures refer to large aggregates 
of atoms. On the other hand, the available data on the mechano-caloric effect 
seem to indicate that the total entropy disappears,* which would mean that we 
cannot separate the entropy of liquid helium into an anomalous and a normal 
part and that even at fairly low temperatures the overwhelming part of the specific 
heat is made up of the excitation of z-particles. More experimental data are clearly 
needed to settle this point. Apart from accurate determinations of the mechano- 
caloric effect, the necessary information can probably be obtained by measuring 
the coefficient of expansion below 1° Kk. The observed decrease in density with 
_ decreasing temperature is an anomalous effect due to “‘ condensation ”’ of z-particles 
and if at the lowest temperatures “‘ lattice vibrations”’ are predominantly*respon- 
sible for the entropy change, the coefficient of expansion should again change its 
sign. , 

The suggested term system removes the difficulty, which has sometimes been 
emphasized in the case of superconductors, that the z-particles appear to be not 
in thermal equilibrium with the rest of the substance. Such a state would of 
course not be stable. ‘The reason why they remain energetically at absolute zero 


* At 1°-5 k., for instance, the observed difference between the entropy of the z-particles and 
the entropy of the whole liquid (4S in formula (2)) is 5x 10-? cal./g., while the total entropy 
has been calculated (Simon, unpublished) as 5-04x10-%cal./g. The extrapolated “lattice 
entropy ” (9=15, Pickard and Simon, unpublished) of 3-8 x 10-2 cal./g. is too big to be accommo- 
dated in the difference between the first two figures. 
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even if the container substance has finite temperatures, is simply due to a lack 
of heat quanta large enough to lift z-particles over the gap. As the temperature 
of the substance is raised, the proportion of large quanta increases until the gap 
becomes insignificant. ‘The shape of the specific heat curves (figure 4) indicates 
that the transition is not uninfluenced by the number of particles in the higher 
State, since it does not follow a Schottky function (figure 64, 1), but indicates a 
cooperative phenomenon like ferromagnetism or the influence of long-range 
order on an order-disorder transformation. ‘The rough term system given in 
figure 7 b does not, of course, allow conclusions to be drawn as to the shape of the 
specific-heat curve, because it is completely empirical and no assumptions are 
made about the density of states above or below the gap. It must be emphasized 
that, since our considerations are based entirely on experimental observations, 
no theoretical explanation for the existence of such gaps in the energy spectrum 
can be given. 
§5. ZERO-POINT DIFFUSION 

So far our account has given no explanation of two outstanding characteristics 
of the z-state, the breakdown of frictionless transport when a critical rate of 
transport is exceeded and the high heat-conduction of liquid helium IJ. Both 
phenomena are, as we shall see, intimately connected. When investigating the 
high heat-conduction of the helium film (Daunt and Mend :lssohn, 1939 a), 
we were able to show that this also was due to a transport phenomenon, the flow 
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Figure 8. The rate of flow surface of liquid helium JI is unaffected by the application of heat to 
the film. This can be explained by the assumption that the z-particles diffuse under their 
zero-point momentum. 


of super-fluid helium atoms towards the place where heat was supplied. We 
then suggested that the heat conduction in the bulk liquid would probably be due 
to a similar though rather more complicated convection process. ‘That this is 
actually the case has since been demonstrated by Kapitza (1940) in a number of 
experiments showing the existence of a counter current of liquid helium in a 
capillary under a temperature gradient. For an investigation of the mechanism 
of the process let us revert to the simple case of the helium film (figure 8). . 
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We observed that the rate at which helium is transferred out of a Dewar 
vessel is not affected when heat is applied to the film outside the vessel. This 
not only demonstrates that the observed high heat-transport through the film is 
entirely due to a transport of helium atoms, but also that this process must evidently 
be isothermal, since a heating of the film at one place would have resulted in 
heat flow into the vessel by conduction. The salient features of the film transfer 
discussed above (§2) show that only z-particles are engaged in the movement and 
the heat supplied to the film must be taken up as heat of excitation which lifts the 
z-particles into the continuum of energy states and evaporates them. Heat is 
of course gained by the helium inside the vessel owing to loss of z-particles, 
but the same gain of heat occurs whether the film is heated or simply allowed to 
drop off. ‘This explanation could be proved by an experiment (Daunt and Men- 
delssohn, 1939b) in which heat was supplied to the inside of a Dewar vessel 
where the concentration of z-particles was thus decreased (figure 9). In this 
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Figure 9. When heat is supplied to the inside of the vessel, the concentration of z-particles is 
decreased and zero-point diffusion takes place along the film from the bath into the vessel. 


case a transfer of z-particles zmto the vessel must take place, resulting, as was 
actually observed, in a rise of the liquid level inside the vessel. ‘The observed 
effect clearly corresponds to the ‘fountain phenomenon”’ discovered by Allen 
and Jones (1938) in the bulk liquid. ‘The comparatively simple conditions of 
our experiment now allow an interpretation of the whole set of effects connected 
with the heat transport in liquid helium II. . 

‘The important point in the experiment shown in figure 8 is that there is no 
“accelerating force”? which drives the z-particles to the place at which heat is 
supplied. It is a striking feature of the experiments that, as was pointed out 
(Daunt and Mendelssohn, 1939 a), a transfer always takes place to those parts 
where the film (i.e. z-particles,) is removed. Since the z-particles do not exchange 
energy with the rest of the substance, the force distributing them evenly over the 
available surface must be their own zero-point energy. The heat transport 
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in a film of liquid helium II is thus simply due to the diffusion of z-particles under 
their zero-point momentum. 

In the bulk liquid, this zero-point diffusion will manifest itself as the tendency 
of the z-particles to take up a statistical distribution within the space occupied 
by the “container substance’’, and a macroscopic movement of z-particles must 
take place whenever this statistical distribution, i.e. the thermal equilibrium in 
the substance as a whole, is disturbed. It may be argued that such a “ diffusion” 
process of a set of particles having zero entropy would be in contradiction to the 
third law. This is, however, not the case since the substance as a whole has finite 
entropy, and the experiments on helium II (Kiirti and Simon, 1938) show that 
at very low temperatures, where the total entropy becomes small, the thermal 
changes associated with the flow phenomenon disappear. It is significant that 
the assemblies exhibiting the z-state—liquid helium and metal electrons—are 
distinguished by their high zero-point energy. The fundamental importance 
of the influence of zero-point energy on the phenomena in helium has been 
emphasized by Simon (1934), who pointed out that its high zero-point energy is 
responsible for the existence of /iguid helium at absolute zero. 

Our explanation of the heat transport in helium II is therefore the opposite 
to the one given by Tisza (1938) on the basis of the hypothesis that liquid helium 
is a Bose-Einstein gas. He attributes mo momentum to the unexcited particles 
and explains the flow of heat as due to an osmotic diffusion of the excited particles 
in the bulk liquid under a temperature gradient. It is difficult to see how such 
a model in which the flow of unexcited particles only appears as a compensating 
displacement current can account for the fountain phenomenon. Quite apart 
from the fact that this model or a similar one proposed by F. London (1939) 
can never lead to zsothermal heat transport, it is in direct contradiction to the 
experiment of figure 8, since there is no bulk liquid for the diffusion of excited 
particles. Such a diffusion of excited particles as visualized by ‘Tisza may of 
course take place as a secondary process in the bulk liquid, and, particularly 
for cases of large heat flow with which the diffusion by zero-point energy cannot 
cope, become the preponderant feature.* This is probably the explanation for 
the peculiar negative fountain effect observed by Allan and Reekie (1939b). 
The inversion point of the fountain observed by these authors will be due to the 
equilibrium between the two rates of diffusion. ‘The striking high heat-conduction 
of the bulk liquid for smail heat flow and in particular its dependence on the 
strength of the heat flow cannot however be explained by ‘Tisza’s model, while it 
agrees well with the process of z-particle diffusion outlined above. 

The conception of zero-point diffusion of z-particles immediately leads to 
an explanation of the critical rate of frictionless transport. ‘Taking again the case 

* The establishment of a temperature difference in helium IJ, therefore, depends on two 
limiting factors, the rate of heat supply and the removal by osmotic and hydrostatic pressure of 
excited particles. The latter is no problem in our experiment (figure 8), where they are simply 
evaporated or drop off, but is important in the bulk liquid. The usual type of heat conduction 
experiments, employing capillaries closed at one or both ends, is thus bound to give results far 
too complex for interpretation, since it is quite impossible to separate the two limiting factors, 


=A more thorough discussion of the heat-conduction process in helium IJ will be given by 
J. G. Daunt and the author later, 
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of the helium film, it is clear that no rise of temperature will take place as long as 
the z-particles which are excited by the heat supply can be replaced by diffusion. 
This rate is limited by the number of z-particles per unit volume (depending on 
temperature) and their zero-point velocity. Once this rate is exceeded, frictionless 
transport breaks down. 

It seems reasonable to assume a similar process of zero-point diffusion for the 
superconductive electrons, since in superconductors, too, we are faced with a 
critical rate of frictionless transport, i.e. the threshold value. When applying 
these considerations to the superconductive electrons we arrive at an interesting 
conclusion. A current set up by zero-point diffusion of electrons must be 
essentially different from a normal current since the momentum of the flowing 
electrons is derived solely from collisions * as in the case of helium atoms. ‘This 
means that the cause of a normal current, the accelerating force of the electric 
field, does not appear at all in the process. This peculiar phenomenon of a 
current in zero electric field is of course exactly what has been observed in super- 
conductors (figure 1 a). According to our conception there should be another 
phenomenon, so far unobserved, by which a super-current is distinguished from 
a normal one. Since it is set up by collisions, the speed of propagation of 
a (heat) impulse should be different, the maximum speed of a signal depending 
on the zero-point velocity of the highest level in A (figure 75). The super- 
current will also be distinguished by having (in a conductor of uniform 
cross-section) uniform charge density. 

Thus zero-point diffusion seems to explain the curious analogy from which our 
considerations originated (1942), i.e. that to the absence of potential differences 
in a superconductor there corresponds the absence of temperature differences in 
the helium film. The conception of zero-point diffusion also removes at once the 
difficulty of understanding the.independence of the transport rate of z-particles 
of the strength of external forces, e.g. the pressure head in helium. It is clear 
that our statement (made in §2) concerning the action of an external force on the 
z-state has to be qualified. Indeed it seems doubtful whether in any of the 
transport phenomena a true particle acceleration occurs. It appears that the 
action of external forces on z-particles is rather that of creating a certain type 
of (space-) order than of accelerating particles. 

To illustrate this point we consider an element (dx) of a long cylinder on the 
surface of which frictionless transport of z-particles can occur. This model 
corresponds to a superconductive wire or to a solid rod covered with the helium 


* When we use here the ordinary kinetic picture of momentum transfer by collision, it is |} 
merely done to distinguish the process of zero-point diffusion from that of acceleration by a macro- _ 
scopic field of force. It has, however, to be remembered that zero-point diffusion is bound to 
differ from an ordinary diffusion process by the fact that while momentum is exchanged, it is not 
dissipated, and it is doubtful whether the terms “ mean free path”’ and “ average velocity ” can 
be used in their accepted meaning. Retaining such a kinetic interpretation of zero-point energy, 
it appears that in systems with high zero-point energy there must be a considerable dispersion of 
velocities. An indication of this seems to have been observed by Ganz (1940) when he found 
that the first indication of a heat impulse travelling through liquid helium was propagated with a 
speed of ~/10* cm./sec., while the speed of the maximum of the impulse was only ~10 cm. |sece 


Incidentally, the secondary maximum observed by him may have been due to Tisza’s diffusion of 
excited particles discussed above. 
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film. If at some place outside (dx), z-particles are removed, zero-point diffusion 
must take place and a macroscopic flow of z-particles in the x-direction will 
occur. ‘This process entails no actual acceleration of z-particles but merely a 
greater mean free path in the x-direction. The flow of z-particles is accompanied 
in the case of electrons by the appearance of a magnetic (but not an electric) 
field and in the case of helium atoms by a moment of inertia. The states of flow 
and rest are thus energetically different but, as the experiments show, of equal 
(zero) entropy. If we make (dx) part of a closed ring, the state of flow, if once 
started, will persist since there are no means of energy dissipation and we arrive 
at a metastable state accompanied in superconductors by a persistent magnetic 
moment and in helium IT by a persistent moment of inertia. The former pheno- 
menon is well known, and it should be possible to detect the latter, particularly 
if these considerations can be extended to a tube filled with the bulk liquid. 

Another way of describing this phenomenon is to say that the momentum 
of frictionless transport is not dissipated because it is zero-point energy. 


§6. SURFACE TRANSPORT 

In our 1942 paper, we have drawn attention to the peculiar tendency of the 
z-particles to move along the geometrical surface of the substance. As a tentative 
explanation it was suggested that this surface transport may be caused by the 
fact that the acting force only penetrates the substance to a certain depth and that 
below this depth no movement of z-particles will take place. It may also be that 
the disturbance of the structure of the “container substance”’ produced by the 
surface facilitates the thermal excitation of z-particles, a process about which 
nothing is known so far. In this case the depth of penetration should be to some 
extent independent of the particular type of disturbance and in some way be a 
measure of the limit of the free path of z-particles. 

In superconductors the phenomenon of surface transport is a well-established 
fact, but in the case of liquid helium the evidence for it is not so strong. In view 
of the similarity of the transport phenomena in the film with those in the bulk 
liquid (Daunt and Mendelssohn, 1939a ; Allen and Misener, 1939) it has been 
suggested that the film may extend below the liquid level and give rise to a transport 
of z-particles along the solid walls in the liquid. Subsequent experiments by 
J. G. Daunt and the author (1939b), Allen and Reekie (1939b), and Kapitza 
(1940) indeed seem to strengthen this view, but it must be pointed out that some 
of the results can be explained equally well by zero-point diffusion of z-particles in 
the bulk liquid. A surface flow of superfluid helium atoms which is analogous 
to the skin effect in superconductors does not however necessitate a higher con- 
centration of z-particles near the walls as is postulated in the model suggested by 
F. London (1939). A decision as to whether and to what extent the frictionless 
transport in liquid helium is confined to the walls of the container has to be left 
to further experiments. 


§7. BOSE-EINSTEIN CONDENSATION 
The conception of the frictionless state of aggregation at which we have 
arrived in the preceding sections was reached merely by conclusions from the 
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existing experimental evidence and does not contain any theoretical hypothesis, 
Even so, a fairly consistent model could be obtained, the chief characteristics of — 
which are those of a quantum liquid condensing in momentum space ; i.e. a set | 
of freely mobile particles of zero thermal energy but an appreciable zero-point | 
energy which is separated from the continuum of thermal states by an energy | 
gap. In 1938 F. London attempted to explain the 4 phenomenon of liquid — 
helium in terms of the condensation of a Bose-Einstein gas, and it is interesting | 
to see how far such a model can be used as a theoretical basis for the interpretation 
of our empirical conception ot the z-state. | 

The chief difficulty is of course that Bose-Einstein statistics is not applicable 
to electrons in metals. However, as has been pointed out by Lindemann (1932), | 
both Bose-Einstein and Fermi-Dirac statistics refer to idealized cases and it is 
quite possible that an intermediate treatment may prove equally applicable to | 
electrons and helium atoms (cf. Daunt and Mendelssohn, 1945). | 

Tt is evident that the high zero-point energy of both liquid helium and the |} 
metal electrons plays a most important part in the behaviour of the z-state whereas 
the (ideal) Bose-Einstein gas has no zero-point energy. ‘This is clearly the reason | 
why Bose-Einstein statistics has failed to provide an explanation for the pressure- | 
independent flow* and the isothermal heat transport in helium films. The | 
process of zero-point diffusion which we have introduced in order to explain 
these phenomena does not exist in the ideal Bose-Einstein gas. Gogate and 
Rai (1944) have recently compared our rate of transfer of superfluid helium 
atoms in films (figure 3 a) with the velocity f of ‘‘non-energetic”’ particles derived 
from Bose-Einstein statistics and have found a fairly good agreement for the 
temperature dependency except for the lower temperatures. In their calculation 
they use the population of the lowest energy state as given by Bose-Einstein 
statistics to be 


wv = nfl (TIT) a eee (3) 


where is the total number of particles per unit volume, 7, the A-point and o 
has the value 3/2. A much closer agreement can however be obtained, as we 
pointed out recently (1945), by using the value 5 for o, and this is the value sub- 
stituted by F. London (1939) in his modified theory of Bose-Einstein condensa- 
tion. It may appear at first as if this remarkable agreement would provide a 
strong support for the application of Bose-Einstein statistics, but actually the 
contrary seems to be the case. ‘The value (c=5) was chosen by London not on 
any theoretical grounds, but merely as the value which happens to fit the observed | 


* Tisza has tried to explain this phenomenon by the accompanying thermal effects, while F. 
London postulates a dependency on the square root of the acting force. Both alternatives have 
been clearly disproved by the experiments (Daunt and Mendelssohn, 1939 a). 

{ In doing this, they tacitly assume the applicability of the quantum-mechanical formula 
v-~h/md for the zero-point velocity of mean free path, where m is the mass of the helium atom | 
and d the mean free path. 'This is very tempting, particularly since adoption of the film thickness 
as value for d leads to a velocity value of the same order as the effective velocity observed in the 
film transport (~10? cm./sec.). The agreement is, however, deceptive, since the zero-point | 
energy of liquid helium calculated by the same formula gives a value 104 times too small. We | 
have therefore adopted the interpretation of the film thickness as a limit (corresponding to the 
lowest levels of A, figure 7b) of free path, given at the end of the first paragraph of §7 of this paper. 
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specific-heat curve of liquid helium.* As was shown in §3, there exists a quite 
general statistical relation which is entirely independent of the particular type of 
statistics employed and which connects the population of the lower energy 
state with the excess specific heat. This relation simply states the fact that the 
population of the lower state determines the change of the degree of order and with 
it the excess specific heat in any thermal excitation of a degree of freedom. The 
case has been thoroughly discussed for order-disorder transformations and for 
molecular rotation (cf. Fowler, 1936). The agreement between the observed 
rate of transfer and equation (3) has therefore nothing to do with Bose-Einstein 
statistics, but merely re-states, and incidentally confirms, the experimental fact 
that the rate of transfer measures the population of the z-state in helium. 

Thus while Bose-Einstein statistics fails to account for the excess specific 
heat of liquid helium and for the characteristic properties of flow and heat transport, 
the question arises whether it has any bearing whatsoever on the z-state as it is 
observed in liquid helium and superconductors. It seems that the position can 
be summarized as follows : the experiments show that both electrons and helium 
atoms can pass into a state in which they are capable of frictionless transport. 
These z-particles are energetically at absolute zero even at finite temperatures— 
they have zero entropy. ‘The existence of particles in the lowest set of quantum 
states does not in itself lead to frictionless transport ; it is necessary that these 
states should be separated from the rest of the energy spectrum by an energy 
gap. ‘That means that a “‘condensation”’ in momentum space must take place if 
the frictionless state of aggregation is to appear, but this need not be a condensation 
into a state of no zero-point energy as in Bose-Einstein statistics. According to 
our considerations (§5) zero-point energy is necessary for the appearance of 
frictionless transport, and in fact no condensation to zero velocity takes place in 
superconductors and probably not in liquid helium either. The undisputed 
merit of Bose-Einstein statistics is that, while it cannot be rigorously applied to 
either superconductors or liquid helium, it is the only statistics so far developed 
which provides for a condensation phenomenon at all ; and by introducing it 
F. London was the first to interpret the A-phenomenon in terms of gas 
degeneracy. 

The problem of liquid helium has been approached more recently by Landau 
(1941) who, instead of introducing a specific theoretical model, considers quite 
generally the quantization of the motion of liquids. He arrives, taking into 
account the experimental results, at two sets of hydrodynamic equations for the 
description of liquid helium. In its present state the theory appears still too 
general to present a rigorous model of liquid helium and it fails to account for 
the essential features of the transport phenomenon, as for instance isothermal 
heat transport, pressure-independent flow and the two components of the fountain 
effect. Its great advantage is that, in contra-distinction to the Bose-Einstein 


* It is interesting that this purely empirical modification, introduced by London to allow for 
the “ unidealty ” of his model, constitutes an approach to our suggested term system (figure 7 b). 
It has been shown that the adoption of (o=5) instead of (o=3/2) seems to produce the 
postulated gap in the energy spectrum of helium between the set of ground states and the 
thermally excited states. 
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model, it does not require modification of an already existing framework but 
seems to be capable of application to the case under consideration by mere 
specialization. 

§8. CONCLUSION 


The picture of the z-state derived from the experimental facts is far too 
consistent to permit us to dismiss the similarity between superconductivity and 
liquid helium II as merely accidental or superficial. Summarizing our conclu- 
sions, we arrive at three salient features of the frictionless state of aggregation : 
(a) The z-particles are energetically at absolute zero. (6) The zero-point energy 
is divided by an energy gap from the thermally excited states. (c) The transport 
phenomena are caused by a diffusion of z-particles under their zero-point 
momentum. 

(a) is fairly well established by direct experiment. ‘The term system sug- 
gested by us to account for (b) is a very rough model and we can hardly expect it 
to yield more than a qualitative explanation, particularly as all evidence for the 
existence of the gap is indirect. (c) is of course closely connected with (a) and 
the fact that frictionless transport is only observed in certain metals is explained 
by (0). (c) is strongly suppported by all the available experimental evidence 
and appears to provide a very satisfactory explanation for the transport phenomena 
in superconductors and liquid helium. ‘The conception of zero-point diffusion 
would seem to be the most important conclusion reached in our considerations. 

The foremost application of these considerations will be to provide a working 
hypothesis for experimental research in a field where so far nearly all important 
progress has been achieved by accidental observations. It is clear that our 
model can be tested by a great variety of experiments, a number of which has been 
indicated in this paper. Apart from the application to experimental work it is 
hoped that the above considerations will also aid the development of a theory 
of superconductivity and liquid helium. ‘The most serious gap in this respect 
is the lack of a kinetic interpretation of a system of freely mobile particles with 
no thermal but finite zero-point energy. 

Finally one may ask whether the z-state is realized anywhere in the world 
except in the laboratory. The high densities in stars and possibly in some of 
the planets may provide suitable conditions for its appearance in macroscopic 
dimensions. Such a high concentration of particles also occurs in the nuclei 
of heavy atoms which have been described by Bohr and Kalckar (1937) as droplets 
of quantum liquid, and Teller and Wheeler (1938) have drawn attention to 
a resemblance between the term system of liquid helium and that of the nucleus. 
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DETECTING TURBULENCE IN LIQUIDS 
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ABSTRACT. The doubly-refractive properties of a weak solution of benzopurpurin | 
were used in a study of the onset of turbulence when the liquid flowed through a long | 
horizontal glass pipe. A beam of polarized light was arranged to traverse a diameter 
of the pipe remote from the ends, the nicols being set to give extinction when the liquid 
was stationary. The emergent light fell on a photo-cell, which was connected through an | 
amplifier to a cathode-ray tube. 

With unidirectional flow, photographs of the trace on the screen showed a steady 
straight line when the velocity was small. Disturbances appeared at a Reynolds number R 
of about 1970, and complete turbulence was established when R attained a value of 2900. 

With forced oscillating flow, the trace at low velocities was a sinusoidal curve, on 
which, at a critical Reynolds number, small superposed ripples were observed. Even | 
when the motion was greatly increased the intricate traces recurred perfectly, but this. 
phase finally broke up into complete turbulence as soon as the amplitude of the surface 
movement in the bottles at the ends of the pipe reached a certain limit. 


§1. INTRODUCTION 


ERTAIN colloidal solutions exhibit the property of streaming doubie 
CC refraction. Most of the quantitative. experiments which have been 
carried out on this phenomenon have been concerned with the problem 

of the shapes of the particles in solution, and they have been performed in the 
well-known concentric cylinder apparatus. More recently, however, this pro-}} 
perty has been used by Alcock and Sadron (1934) to provide a means of measuring} 
the velocity distribution in various kinds of laminar flow. The properties of] 
many liquids, suitable for the purpose, have been examined by Weller, Middlehurst} 
and Steiner (1942), who gave numerous references to previous work. Andrade 
and Lewis (1926), who used rotating cylinders, remarked that the onset of tur- 
bulence was favoured by the presence of colloidal particles, but apart from al 
brief note by Hauser and Dewey (1939), no trace has been found of any attemp 
to use the method deliberately for detecting turbulence. However, the under-} 
lying principle appeared promising for the purpose and worth detailed attention.] 
Accordingly both unidirectional and forced oscillating flow in a circular pipe 
were examined by this means at the Engineering Laboratory, Oxford. Many, 
investigators * agree that, in the former, laminar flow breaks down at a Reynolds} 
number of about 2100, calculated with the pipe diameter and with the mean 
velocity of flow. On the other hand, no previous information is forthcoming: 
concerning forced oscillating flow. 


* Cf, Goldstein (1938). 
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§2. THE WORKING LIQUID 


The best known colloidal solutions possessing the required properties are 
vanadium pentoxide (V,O;), the red dye benzopurpurin 4B and certain clay 
suspensions. ‘The first two contain filamentous particles ; when stationary, 
the filaments are oriented at random owing to thermal agitation, but when in 
slow motion they cause the liquid to exhibit doubly refractive properties if viewed 
at right angles to the direction of flow. It appears from a review by Edsall 
(1942) of the present state of knowledge that there is a general measure of agree- 
ment amongst physical chemists concerning the nature of this effect. In steady 
laminar flow, the filaments as they are carried with the stream rotate about axes 
perpendicular to the direction of streaming, but they do so with an angular 
velocity which is a minimum when they lie along the stream. Thus at any instant 
there will be more particles in this position than in any other. This explanation 
also shows why these solutions, unless highly diluted, exhibit non-Newtonian 
effects, such as the distortion, examined by Lawrence (1935), of the parabolic 
velocity distribution in laminar unidirectional flow in a pipe. For it is clear 
that the filaments, as they whirl along, must be comparatively far apart if no 
interaction is to occur. 

Thus for the present purpose it was necessary to select a solution sufficiently 
weak to show no appreciable anomalous properties, yet strong enough to exhibit 
measurable optical erfects. After some preliminary tests, 0-25°,, benzopurpurin 
solution was chosen because it combines ease of preparation in quantity with a 
comparatively powerful double-refraction effect. ‘The solution was prepared 
2 |. at a time ; 0-5 gm. of the powder was ground up in a mortar with a little 
distilled water, and the volume was made up to 0:5 1. After a slight warming 
the powder was completely dissolved, giving a clear red solution which was then 
diluted to21. The alternative procedure of adding the powder to the full volume 
of water was found to be ineffective ; the powder then refused to dissolve even 
on boiling. A solution of 80 gm. of potassium sulphate in 900 cc. of distilled 
water was then made, and 200 cc. of this was added to each 2 |. of benzopurpurin 
to form the working liquid. The salt solution causes the dye molecules to aggre- 
gate to particles of a size suitable for the double-refraction phenomenon, but 
larger quantities throw out the dye as a flocculent precipitate. ‘The aggregation 
is caused by the electrical effects of the monovalent kation ; the presence of 
kations of higher valency, derived, for example, by corrosion, must be guarded 
against in view of their powerful coagulating effect. Certain parts of the 
apparatus had to be made of metal, therefore the sulphate was used in preference 
to the more corrosive chloride. 

One disadvantage of benzopurpurin as a working liquid is the change of its 
properties with time. For a few hours after mixing it was found to be inactive. 
On the first day after mixing its optical power was sufficiently strong for the 
experiments to proceed, but this increased too much during the day’s work 
for the results to be entirely satisfactory. To avoid this difficulty, observations 
were then confined to the third and to the seventh days, over which the rate of 
levelopment was slow. The growth of the particles depends upon the tem- 
yerature, but this could not be kept constant ; hence the optical properties of 
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the liquid on the two test days was not the same in successive sets of experiments. 
The liquid remained active for periods longer than a week, but it gradually 
became cloudy and therefore unsuitable. 

When it was desired to interpret the observations in terms of Reynolds 
numbers, considerable difficulty was experienced because the kinematic viscosity 
of the working fluid was not free from uncertainty. Mark (1940), who described 
a benzopurpurin solution as “elastic”, explained that its viscosity is maximum 
at very low shearing stresses, diminishing in an exponential manner to a steady 
Newtonian value as the shearing stress is increased. ‘These effects were dis- 
covered in the working fluid, although it was so dilute. Experiments were made 
with a viscometer of the Couette type, in-which the diameter of the outer rotating 
cylindrical surface was 12 in. and that of the inner cylinder (hung on a torsion 
wire) was 2in. Ata mean velocity gradient of 3 sec. the viscosity of the three-day- 
old fluid was found to be 36% in excess of that of water, while at 9 sec.+, the 
limit of the apparatus, the corresponding value fell to 21%. At a very high 
velocity gradient, estimated at 800 sec.—!, observations with an Ostwald (capillary) 
viscometer gave the excess of kinematic viscosity as only 3%. The experiments 
to be described lay in a comparatively small intermediate range, and the onset 
of turbulence in unidirectional flow was at a mean velocity gradient of about 
17} sec... ‘lhe density of the working fluid was }°% greater than that of water. 
On the assumption that the excess of kinematic viscosity diminished exponentially 
with the mean velocity gradient, logarithmic plotting of the above results showed 
that at 174 sec.-! the excess was about 9%, and, as no sensible change could be 
detected after seven days, this figure has been used in all the calculations. 

It should be borne in mind that the viscosity of these anomalous solutions 
is known to depend partly on the previous mechanical history of the solution ; 
vigorous treatment breaks up temporary linkages of the particles and for a time 
reduces the viscosity. Now in the unidirectional flow experiments the working 
liquid passed continuously through a centrifugal pump, and the oscillating 
flow apparatus was operated as violently as possible for a few minutes at the 
commencement of each series of tests. ‘Therefore viscosities deduced from 
measurements made under the quiet conditions in the Couette viscometer were 
probably somewhat higher than those obtaining in the experiments. Support 
for this view is provided by the unidirectional flow tests described in §4. There 
it was found that the critical velocity occurred at a Reynolds number of about 
1970, which is slightly lower than that established by other methods. 


§3;- THE OPTICAL APPARATUS 


The horizontal glass pipe, through which the liquid passed, was 15 ft. long. 
It was made up of three 5-ft. lengths accurately fused together, the mean internal 
diameter of the central length being 0-94, in. At the central section, distant 
95 diameters from the ends, the flow was examined by means of a photo-electric 
cell. ‘The source of light was a 125 w. tungsten arc lamp (pointolite), the hori- 
zontal beam from which passed successively through a convex lens (to make 
the beam parallel), a cell containing a solution of ferrous and copper sulphates 
(to absorb injurious heat rays), an adjustable stop set at about 7/sp in. diameter, 
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the first nicol or polarizer, a diameter of the pipe, and finally the second nicol 
or analyser. The axes of the nicols were fixed at 90° to each other, both being 
at 45° to the axis of the pipe. In order to preserve the approximately parallel 
nature of the beam,a rectangular glass cell full of distilled water was fitted round 
the section of the pipe which was under observation. 

When the liquid was stationary, with the dye filaments oriented at random, 
the only light which emerged was that due to lack of parallelism in the beam, 
causing unequal reflection losses of the components of the plane-polarized light 
along and across the tube, with consequent imperfect extinction by the nicols. 
But, on putting the liquid into motion, the double-refraction effect caused the 
emergent beam to increase in strength. It fell on a photo-electric electron 
multiplier, consisting of a caesium photocell whose primary emission was mag- 
nified about 2:5 x 104 by seven stages of secondary photo-emission controlled 
by crossed electromagnetic and electrostatic fields. It is free from inertia effects 
in rapidly varying light. ‘The output of the multiplier was led through a single 
stage of valve amplification, having a voltage gain of about 20, to the Y plates of a 
gas-focused cathode ray tube, the screen of which had a long after-glow. ‘Thus, 
even at the very low velocities of sweep which were used, it was possible to compare 
the traces of successive cycles. ‘lhe time base was of orthodox type, a condenser 
being charged through a pentode valve and discharged by contacts mechanically 
operated by the piston rod of the reciprocating pump described in §5. In 
order to secure the necessary stability, all potentials were supplied by dry batteries 
and accumulators. With a Leica camera, 56 photographs of unidirectional 
flow were obtained and 123 of oscillating flow ; those reproduced in figures 2, 
5 and 6 were taken at approximately half gain, which was found to be the most 
suitable ratio because at full gain, in spite of many precautions, the trace was 
disturbed by stray fields. ‘The zero level of the trace was indicated by projecting 
the images of spots or arrows on to the screen. 

In order that a satisfactory comparison might be made between the two types 
of flow, no adjustment was made to the glass tube or to the optical apparatus in 
the course of the work detailed in the following pages. Owing to the difficulty 
of obtaining a circulating pump, the experiments on oscillating flow were performed 
first. Here, however, the work will be described in the reverse but logical order. 


§4. THE BREAKDOWN OF UNIDIRECTIONAL LAMINAR FLOW 


As shown diagrammatically in figure 1, the inlet end of the glass pipe was 
joined by hose and a short length of 1-in. stainless-steel pipe to an aspirator 
bottle A. The outlet end was similarly connected to glass tubing leading to the 
large bottle B ; this served as a reservoir from which the pump C drew its supply. 
The pump was of.the centrifugal type, driven by an electric motor; its interior 
was chromium-plated and, to give additional protection, the plating was coated 
with bakelite varnish. ‘The delivery from the pump passed either to the bottle A 
through glass tubing and an 11-ft. length of 1-in. stainless-steel pipe or back to 
the bottle B through a by-pass. Thus the flow through the glass pipe was tur- 
bulent at entry ; it was effected by gravity only, and any pulsations developed in 
the pump could not reach the section under observation. ‘The flow was con- 
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trolled by means of pinch-cocks on the by-pass and near the outlet end of the glass 
pipe, and also by altering the speed of the pump. No difficulty was experienced 
in maintaining constant conditions. The apparatus required a charge of about 
8 1. of working fluid. 

To measure the flow, a stainless-steel orifice-plate was inserted at D, and open 
glass stand-pipes fixed to a vertical scale were connected to both sides of it. In 
view of the limited supply of working fluid, it was thought preferable to use tap- 
water when determining the coefficients of discharge of the orifice at various 
rates of flow. In order that the observers should not be biased in any way, this 
calibration was not done until all the photographs had been obtained. 

Observation revealed, and the camera recorded (figure 2), that at low velocities 
the trace on the screen was a horizontal straight line, which moved upwards 
further from the zero level when the velocity was increased. ‘The zero with the 


Figure 1. Elevation of the unidirectional flow apparatus. 


tungsten arc on and the liquid stationary is indicated on the photographs by 
two spots ; when the light was turned off, the trace descended only by an amount 
approximately equal to its own breadth. ‘The sweep was set to operate at the 
constant speed of 20c./min. At higher velocities of flow, momentary oscillations 
occurred, evidently due to the “‘ flashes” described by Reynolds. As the velocity 
was still further increased, the disturbances became more numerous and the 
frequency of the ripples went up. Finally a velocity was obtained when the 
trace was always in a state of agitation. 

With the three-day-old liquid, photographs were secured of complete steadi- 
ness at a Reynolds number R=1970, and of disturbances at R=2270. But a 
closer “bracket” was obtained after seven days, therefore a selection of photo- 
graphs from this test is reproduced in figure 2, where (i) and (ii) show the steady 
straight lines at Reynolds numbers of 750 and 1900. At R=1970 there were 
two “flashes” in five minutes, but attempts to photograph them were unsuccessful ; 
however, a picture (iii) taken during the intervening periods shows that the steady 
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(11) R=1900 


(iii) R=1970. (iv) R=2030. 


(v) R=2620. (vi) R=2910. 


Figure 2. Breakdown of laminar unidirectional flow. 
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line was slightly bowed. At R=2030 the disturbances were marked, a period 
of rippling and its termination being shown on (iv). At higher velocities the 
quiet periods were shorter and the rippling more violent, although even at R= 
2620 (shown in v) the bursts of turbulence did not often last longer than two 
cycles. In this phase the trace in an inactive period was seldom perfectly hori- 
zontal, and the approach of a turbulent spell could often be detected by a gentle 
fallof the line. At R=2910(vi) really quiet intervals never occured, and turbu- 
lence was completely established, while at still higher velocities the amplitude 
and frequency of the ripples further increased. 

These results show at once that the apparatus gives much the same value of the 
critical Reynolds number as that obtained by earlier methods, but further informa- 
tion may be obtained from the photographs. Figure 3 gives the distance of the 


zy 1000 R 2000 3000 


Figure 3. Deflexion of trace in unidirectional flow. 
o After 3 days. e After 7 days. 


straight traces from the zero with no flow plotted against R. Little reliance 
can be placed on points beyond the critical value because, as explained above, 
the line in a quiet period was not often truly horizontal, and these measurements 
are therefore dubious. Confining attention to laminar flow, we see that in the 
three-day tests the deflexion was linearly dependent on R, and this suggests 
little orientation of the particles at the very low tangential stresses set up 1n such 
a wide pipe. Had the orientation been complete, the curve would have become 
horizontal ; in the seven-day tests it did in fact tend to do so, but the onset of 
turbulence cut off the finish of the process. These ideas are confirmed by the 
observation that, when the flow was stopped, the deflexion of the trace did not 
immediately fall to zero, from which it follows that the particles cannot have been 
very small. ‘Thus the low degree of orientation must be attributed either to lack 
of anisometry in the particles, which seems improbable in view of their known 
filamentous shape, or to the small tangential stresses. 

The investigation having for its main object the detection of turbulence, no 
attempt was made to analyse the turbulent vibrations, but it may be remarked 
that the laminar flow evidently broke up into vortices not much smaller than the 
cross-section of the light beam used. Only if the vortices were considerably 
less than this would the double refraction be reduced to zero, and this was not 


observed. 
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§5. THE BREAKDOWN OF OSCILLATING LAMINAR FLOW 


It is desirable first to examine theoretically the nature of forced oscillating 
laminar flow ina pipe—a problem which has been attacked by Southwell and 
his collaborators, e.g. Christopherson, Gemant, Hogg and Southwell (1938). 
They considered a U-tube (composed of a horizontal pipe connected at 
each end to a bottle) in which the motion was maintained by a pulsating 
air pressure. Being concerned principally with establishing the basis of their 
method of viscometry, they concentrated their attention on the relation 
between the amplitudes of the applied pressure and of the surface movement in 
the bottles. Here, however, we are interested in the distribution of horizontal 
displacements over a cross-section of the pipe which results from a measured 
amplitude of the surface movement in the bottles. Accordingly Southwell’s 
analysis has been extended. ‘The details are relegated to the Appendix, but the 
results are shown in figure 4. ‘This is a vector diagram drawn for the motion 
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Figure 4. Theoretical distribution of displacement in laminar oscillating tiow. 


in the 15-ft. pipe of water at 15° c. at 8 and 20 c./min., these frequencies being 
the slowest and fastest used in the experiments. The direction of the line Ox 
represents the amplitude in the bottles. The curve is graduated in terms of the 
ratio of the radius r under consideration to the radius a of the pipe, and the 
magnitude and phase angle of the horizontal displacement are given by the straight 
line drawn from O to the curve. The distribution of these displacements has a 
kind of inverse analogy to the well-known skin effect which occurs when high- 
frequency alternating current passes through a solid circular conductor. It will 
be seen that at 8 c./min. the liquid within a radius of about 0-44 oscillates with 
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little relative movement. At 20 c./min. this radius is increased to about 0-6 a : 
it is only at frequencies beyond the scope of the experiments that almost the 
whole of the liquid oscillates virtually as a solid plug with large shearing stresses 
close to the wall. It should be noted that the ratio of the horizontal displacements 
in the pipe to the amplitude in the bottles is independent of the latter when the 
other variables are unchanged. 

In the first set of experiments the glass pipe was connected at each end to an 
aspirator bottle of diameter 5-6 in. In one bottle, off which the top was cut, the 
movement was measured with a stainless-steel hook gauge, surface disturbances 
being eliminated by perforated bakelite baffles. The other bottle was joined 
by rubber tubing to a single-cylinder reciprocating pump of diameter 3 in. 
and stroke 13 in., from which the valves were removed. The piston was actuated 
by a Scotch crank driven through reduction gearing by an electric motor, so that 
its motion was truly simple harmonic. It was found that if, on starting, the rubber 
tubing was connected when the piston was at the middle of its stroke, no creep in 
the excursions of the liquid took place. The amplitude of the motion of the 
liquid was controlled with four bottles (serving as capacities), joined singly to 
the connexion between the pump and the aspirator bottle; they could be turned 
on or off at will, and the effective volume of one could be finely adjusted by 
partially filling it with water. ‘To satisfy the requirement that the pressure 
imposed on the working liquid should vary in a sinusoidal manner, the pressure 
changes in the air were kept small and nearly isothermal by running the pump 
only at very low speeds. ‘lhe apparatus required about 5 1. of working liquid. 

‘The method of test was to take a series of photographs at constant pump speed, 
the amplitude of the movement in the aspirator bottles being varied by means of 
the capacities. ‘The surprising fact at once emerged that over a considerable 
range the trace at each amplitude, although considerably agitated, recurred 
exactly for an indefinite length of time. ‘The break-up of this phase at 13 c./min. 
is illustrated by the series shown in figure 5, in which are indicated the travels 
measured in the aspirator bottle and the Reynolds numbers. ‘The latter will 
now be formed with the maximum velocity of the motion, with the simplifying 
assumption that the velocity was uniform over the cross-section. The zero 
at no flow is here indicated by the tips of the arrows. ‘The first photograph 
shows recurring flow at a travel of 0-55 cm. ; the ends of the stroke correspond 
to the undisturbed portions of the trace, which, as in the unidirectional experi- 
ments, did not instantly descend to zero when the velocity became very small. 
At 0-80 cm., (ii), the trace still recurred. ‘The next photograph records two con- 
secutive cycles, the only difference which can be detected between them being 
the short-circuit of the prominent inverted W on the right. Again, the ends of 
the stroke may be distinguished by the quieter motion there, which, being slower, 
caused a thickening of the trace. Finally, at 1-61 cm., (iv), the diagrams did not 
recur, although the differences between them was not marked, but at still-greater 
travels the motion of the trace was wildly irregular. ‘These results indicate 
that under certain conditions a recurring system of vortices was formed, which 
was destroyed when the motion became too violent. ‘The available information 
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concerning the region of breakdown is summarized in table 1, in which the 
results of the one-day test (mentioned in § 2) have been included. 


Table 1. The breakdown of recurring turbulent flow 


Age of fluid — c./min. ‘Travel (cm.) R 


One day 20 0):97 6360 Recurring 
a 4 1-30 8480 Almost recurring 
re i 1-41 9220 Not recurring 
= 164 1:24 6840 Recurring 
f ak 1:52 8360 Not recurring 
a 13 1:03 4610 Recurring 
; % 139 5880 Almost recurring 
re " 1°51 6750 Not recurring 
Three days 20 0-98 6500 Recurring 
7. in 1-32 8800 Not recurring 
; 13 1:08 4840 Recurring 
3 , 1:39 6190 Not recurring 
Seven days 20 0-95 5810 Recurring 
Se rr ils 8050 Not recurring 
- 13 0°80 3480 Recurring. 
yp 3 1-10 4790 Almost recurring 
9 5 125 5460 Almost recurring 
96 y iho@il 6990 ~=Not recurring 


No correlation can be discovered among the values of R, but there are indications 
that it was the travel which was the determining factor. Allowing for the difh- 
culty of estimating when recurrence was perfect, it may be suggested that a travel 
of about 1-25 cm. was the criterion of breakdown. 

Attention was then concentrated on the first signs of turbulence. With a fresh 
charge of liquid, the trace at very slow flows was seen to be a smooth sinusoidal 
curve indicating laminar conditions. As the amplitude was increased, small 
ripples appeared at first, but later the harmonic nature of the picture was broken 
up into the form indicated in figure 5 (i). ‘These experiments were found difficult 
and unsatisfactory because under some conditions the travel did not remain 
steady. Consequently suspicion fell on the apparatus used to maintain the 
motion. It has been pointed out by den Hartog (1940) that, if a forced oscillation 
is maintained with a spring of varying stiffness as coupling, the resulting motion 
may be unstable. In the experiments, the damping was considerable owing 
to the destruction of kinetic energy when the liquid entered an aspirator bottle ; 
free oscillations were in fact nearly dead-beat. Thus the power required to 
drive the apparatus, and consequently the compression and varying elasticity 
of the intermediate air, seemed unnecessarily large. Moreover, at 20 c./min. 
the travel required for laminar flow was too small to be measurable with reasonable 
percentage accuracy ; and when the frequency was halved in order that a wide 
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(1) Travel 0-55 cm., Max. R=2400 (ii) Travel 0:80 cm., Max. R=3480. 


(iii) Travel 1:25 cm., Max. R=5460. (iv) Travel 1-61 cm., Max. R=6990. 


Figure 5. Breakdown of recurring oscillating flow. 
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(i) Travel 4:39 cm., Max. R=730. (ii) Travel 5:83 cm., Max. R= 980 


(iii) Travel 6:25 cm., Max. R=1050. (iv) Travel 8-22 cm., Max. R=1360. 


Figure 6. Breakdown of laminar oscillating flow. 
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range of conditions might be investigated, it approached so closely to the natural 
period that sufficiently small travels were unobtainable. Accordingly the 
aspirator bottles were replaced by vertical glass tubes, 19 in. long and only 13 in. 
in diameter, the ‘ower portions of which were bent through an easy right angle. 
A carboy was substituted for the four capacities, and the resistance to the motion 
of the air was reduced by increasing the diameter of the connecting tubing from 
4 tozin. The travel was measured with a horizontal telescope mounted on a 
cathetometer. 

As expected, these modifications greatly reduced the damping, and tests 
were made at frequencies both above and below the natural period, but the 
difficulties mentioned above did not disappear. A test after three days at 20 c./min. 
showed, after a considerable wait, a smooth sinusoidal curve at R=890. At the 
next stage, a recurring trace with slight superposed ripples was obtained after 
allowing it a long time to settle, but the travel never became absolutely steady : 
its mean was 2:51 cm. (corresponding to which R is 1050), but the excursions 
varied by about } mm. at each end. ‘These irregularities in the travel continued 
up to R= 1320, accompanied by rippling of increasing violence, but at and beyond 
R=1480 no trouble was experienced in quickly obtaining both recurring traces 
and steady travels. Other observations made it clear that it was the Reynolds — 
number, not the travel, which was the criterion of the onset of turbulence and of 
irregular travel. This conclusion is supported by the illustrations displayed in 
figure 6, which were taken with fluid seven days old at a speed of only 84 c./min. 
Although the sweep unfortunately gave trouble at this point, (i) shows that at 
R=730 the trace was smooth. ‘The displacement of the trace being greater than 
in the three-day tests described above, at K=800 and 920 it was just possible to 
detect faint ripples, but they are too small to be seen on the necessarily small 
scale to which the photographs are reproduced here. At R= 980, (ii), the rippling 
was more marked ; at R=1050, (iii), however, the difficulty of obtaining steady 
travel again commenced, although the trace recurred. he irregularities were 
at their height at R=1240: here the trace apparently settled down quickly, 
but then a spasm occurred and the trace slowly changed its form, not settling 
down to its original shape until 20 min. had elapsed. ‘The recurring diagram, 
(iv), at R= 1360 was taken after 15 min., the last 5 min. being steady, but at and 
beyond R=1610 no further trouble was encountered. 

The results suggest a critical value of R, not very sharply marked, at about 
800. The existence of a zone of irregular travel at slightly greater values of R 
may be attributed, not to the drive, but to the formation of small eddies, prin- 
cipally at the ends of the pipe. ‘These would erratically increase the small 
viscous resistance to motion, and, when sufficiently large, would cause a noticeable 
variation in the travel. At greater amplitudes they would become more numerous 
and possess an almost constant mean value, and under these conditions a steady 
travel would again be obtained. It thus appears that the onset of turbulence 
can be detected, although in a less sensitive manner, simply by examining the 
uniformity of the travel, 
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§6. CONCLUSIONS 


(i) Considerable difficulty was encountered in securing a satisfactory com- 
promise between the conflicting requirements of the various parts of the apparatus. — 
A wide pipe and a narrow beam of light were employed so that the lens effect of | 
the former was small. Thus the emergent beam was sensibly parallel, and it 
was possible to obtain good extinction with the analyser. The disadvantage of 
the wide pipe was a very low critical velocity, occurring at such small shearing 
stresses that a slight anomaly in the viscosity of the working liquid, with its 
attendant distortion of the velocity distribution, was present. If, to avoid these 
effects, the working liquid had been a weaker solution, even more sensitive elec- 
trical apparatus would have been necessary. 


(ii) With unidirectional laminar flow, the trace was a horizontal straight 
line. Ata Reynolds number of 1970, brief disturbances appeared, which became 
more numerous and violent as the velocity was increased. ‘Thus the indications 
of the apparatus were in fair accordance with the results obtained by other 
methods. ‘Completely developed turbulence was shown to set in at a Reynolds 
number of about 2900. 


(iii) With oscillating laminar flow, the trace was a smooth sinusoidal curve. 
At a certain Reynolds number (derived in an approximate manner) recurring 
ripples appeared, which suggested the formation of a steady system of vortices. 
Over a small range of slightly higher Reynolds numbers the motion of the liquid 
became slightly erratic owing it is thought, to the formation of eddies which 
irregularly increased the viscous resistance. Under much more vigorous con- 
ditions of flow the recurring traces broke up, indicating complete turbulence, 
and this point was found to depend upon the travel of the liquid. 
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APPEN DIX 
Oscillatory horizontal displacements in the pipe 


An expression is required giving the distribution of horizontal displacements 


over the cross-section of the pipe in terms of the amplitude of the surface move- 
ment in the bottles, 
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We adopt the notation used by Southwell and his collaborators (loc.-cit.), 
and refer to their results by inserting the numbering of their equations at the left- 
hand margin of the page. Let p be the externally applied pressure, x the displace- 
ment in the bottles, v and y the horizontal velocity and displacement in the pipe 
at radius r. ‘Che uniform cross-section of the bottles is denoted by A, a is the 
radius of the pipe, w and p are the viscosity and density of the liquid, 2h is the 
length of the pipe, and g is the acceleration due to gravity. 

We assume that 


(? as Pe, ) 
4 x= Xen, | 
@) i: v= Vent > cig Fe (1) 
y= Yeint, J 


where P, X, V and Y are complex, and V and Y are functions of r only. Then 
it was proved that 


ual eee J o(kr) _ 
(8) V = Re @ an) E Taal Niecerstctcreis (11) 
OREM ENE 
5) { where O= Toe | rVar, 
and k?2 = — imp/p. 
hu 
De Dut X =—Q, 
(9) Bu a Q 
res = J,(Ra) oe 8, fe: 
Beene i {! ka Jo(ka) g 2phn®’ 
,_ 7a" 
(7) where = a | 
Hence, on integrating (i1), we obtain, after some reduction, 
yg Z J.(Ra)\= (1 Ga) ee 
2 =<1— = settee MU ARS: 
X¢g i! ha J (ka) J o(ka) a 


where F£ is a function of r only. It is evidently zero when a steady state of oscil- 
lation has been attained. Now k? is a pure imaginary, hence we introduce a 
quantity z given by 
(12) gia — kitgt, 
Then (iii) may be written 
Yee elie He (1 ev). a. (wv) 
izv/t J ((izv/2) 


J (iz v/2) 
Here z and w are real, and they are given by 


$ + 
z= a(2), war(*) 
we be 
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The evaluation of the right-hand side of (iv) may be performed in two parts. 
The first term reduces to 
( C(z) -1D(z), " 

| A(z)=1 HO 


where 2)= So zs) = | 
a0) 4 CON Gay BaP OO" AGT Be! > 
_ x We) _ & XZ) 


The modulus of this vector is [{C(z)}?+ {D(z)}*]? and its phase angle @ is given 
by tan @= —D(z)/C(z). Numerical values of these expressions can be obtained | 
with the aid of the tables given by Russell (1914). The second term in (iv), | 
which gives the displacement distribution over the cross-section, may be expe- | 
ditiously dealt with by writing it in the form 


— Ml) 51:40,0)—.+ 
M, (2) 
and making use of McLachlan’s tables (1934). 

The diagrams of figure 4 can then be constructed. ‘Taking OX to represent 
the direction of the vector X, we set off OP of length [{C(z)}*+ {D(z)}?]* at the 
phase angle 6. For various values of r/a, rays from P are drawn of length} 
[{C(z)}2 + {D(z)}]? My(@)/M(z), the angles between them and OP being] 
{8)(w) —O(z)}. Lastly, the curve joining the ends of the rays is drawn in. 
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ABSTRACT. From a pair of axial magnification formulae, exact addition equations are 
derived in this paper for the axial spherical aberration and curvature of field for central 
and sagittal rays produced by any number of centred spherical surfaces whatever their 
individual separations. 

For a ray of initial semi-aperture U, and final semi-aperture Ug’ passing through media 


of refractive indices m, mg... . mq, the condition for zero axial spherical aberration is 
p=1 
XL Npupcy (sin Up + sin Ip— sin Ip’ — sin Up’) = 0, 
p=a 


where cp is the distance of the pth intermediate paraxial object point from the ‘centre of 
curvature of the pth surface. 

For zero Petzval curvature and an initial semi-angle of field 0, the condition is first 
eiven in the usual approximate Petzval form with the addition of a corrective term 


P=1 ny’ —n 0 Ny x 1 2 
ee ear a 2— mp a tan= Op =*() 
p=q Up Nprp Np D 2 


jute 
A more simple alternative form, however, is X 8p tan $ 4)—O0, where dp is the deviation 
p=4 


of the central paraxial ray at the pth surface. 
An addition equation is also given for the sagittal astigmatism. 


§1. INTRODUCTION 


xAcT values for the aberrations of the image formed by compound optical 

systems of centred spherical refracting surfaces are usually obtained by 

trigonometrical ray-tracing, and the removal of the aberrations, although 
upplemented by well-known approximate formulae in many cases, thus depends 
n some degree on a trial-and-error process, in which values for the unknown 
adii of curvature are first assumed and then modified or adjusted until the 
elected aberrations are removed. 

The object of this paper is to show that so far as the axial spherical aberration, 
he curvature of field for central rays (Petzval curvature), and to a lesser degree 
he sagittal astigmatism, are concerned, exact addition formulae can be derived, 
ontaining a minimum number of variables, which, when equated to zero, give 
1e direct conditions for the absence of these aberrations. 

It is hoped that the publication of these formulae may, in combination with 
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the well-known use of the sine condition for the pre-determination of coma, 
indicate useful means of estimating in advance the magnitude of the several 

aberrations in any given system, and at the same time encourage further develop- | 
ment along these lines. 


§2. THE AXIAL MAGNIFICATION FORMULAE 


The formulae developed in this paper are based upon a pair of axial magnifi- | 
cation formulae connecting (a) the successive intercepts A/ and Al’ cut off 
respectively before and after refraction by a pair of paraxial rays, and (4) the’ 
corresponding intercepts AL and AL’, when one of the pair is paraxial and the 
other of finite semi-aperture U. 

For infinitesimal lengths d/ and dl’ it is well known that the relation for. 
paraxial rays is | 


dl’ =~ m2dl, 
n 


but for intercepts of finite length, this equation is an approximation only. 
The exact equation of transfer for a pair of paraxial rays is 


/ 


Al’ = —mmAl ors ¢u'uy Al =n AL |) ee ee (1) I 


where m and m, are the linear magnifications at the respective extremities of the } 
length Al’, wu and up are the angles between the respective paraxial rays and the} 
axis before refraction, and w’ and uw’, are the corresponding angles after refraction. || 
This equation may be deduced from the familiar paraxial formulae 


If =—sg'=—80o8 and fig=m f'/gs=n'm,/n, 
remembering that Al’=g’ —g,’, where the distances g, 9, g’, 2) are measured| 
from the principal foci f and f’ respectively. 


It is sometimes convenient to substitute m for mp, by putting g+Al for go, so 
that 


Al 
m _ gt Seb 2s nee 
Mp g fi 
whence 
1 Areas 
My m : fae 
similarly 
cei 
i tty. fo. (2) 


Substituting the value of m, above in equation (1), 


Al = = mel (1 2 “~—), 


i 


an equation which will be useful later on in the paper. | 

Next, let one of the rays be of finite aperture U cutting the optical axis in 
and B’ before and after refraction at distances C and C’ from the centre of curva} 
ture of the refracting surface (see figure 1). Let b and b ‘denote the correspondin 
intersections of the paraxial ray at distances c and c’ from the centre, Final 
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let b) represent the position of the paraxial image of B, and c,’ its distance from 
the centre of curvature. 


B 
Figure 1. 
OBE C—c. Bb—Ay Al: 
OBZ Cs By oy VA EY, 
Ob—« by’b’=4I’, 
Ob’=c’, paraxial image of b. B’b’= A=longitudinal spherical 
Obo’ =p’, paraxial image of B. aberration. 


Then, using AL’ and AL in place of Al’ and Al in order to indicate that one 
of the image-forming rays is now of finite aperture, and denoting the magnifica- 
tions at b’, by’ and B’ by m, my and M respectively, we have, from equation (1), 


ING ~ mm, AL +A, 


where A is the longitudinal spherical aberration of the image of the point B for 
semi-aperture U. 
Substituting cy’/C for my in the preceding equation, we have 


INES OG 
TG 


ee a 
n (€ 
n' 


nm AL Se ON 
= mMAL+A(1 +P = me), since 7G = Mo: 


n 


id 


- 5 n 
Again, since — mmAL=ATl', 
n 


nt AL _ cy Al 
Ne? ly ee 
Hence ; 
Aid = Ren hae iy elie ee Ege (3) 
n Co n'u’ sin U Co 


In a system of centred spherical refracting surfaces, the AL’ of one surface 
becomes the AL of the next, and the U’ of the one surface the U of the next. 
Equation (3) may therefore be extended to 


p=1 Cc / 
(Shae , (pins u eae Has , Dp 
nu sin UAL =nyu, sin U,AL, + poe u,' sin U, AD fides (4) 
7% ‘Pp 
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If in equation (4) a substitution can be made for the A, in terms of the | 
known elements of the pth surface, we have available an addition equation for | : 
spherical aberration. By a happy accident, as will be seen in the next section, || 
we can not only substitute a simple expression for A,, but also get rid of all the | 
troublesome quantities of the type cj,. i 


§3. SPHERICAL ABERRATION OF A RAY OF FINITE APERTURE JU, |i 
PRODUCED BY A SINGLE SPHERICAL SURFACE, SEPARATED BY ]if 
MEDIA OF REFRACTIVE INDICES n AND wn’ i 
In the present section, departing from accepted procedure, the radii of | 

curvature of the refracting surfaces are treated as unknown quantities to be | 

subsequently determined, and the formula for spherical aberration is based, in | 
the first instance, on the semi-apertures U and U’ before and after refractions. | 

‘These are selected as known quantities and may be said to determine the shape 

of the system. || 
It only remains to fix a single linear constant (conveniently represented by || 

paraxial c’) to determine the scale, after which all other linear quantities are | 
uniquely determined, including the curvature of the surface and the spherical |] 
aberration. | 
By definition we have 
A C-c .nsnU mu ni fae —sm UU 
nu’ ral sin U’ i. 


if sin U=au. 
If u be taken so that a=1, the preceding equation becomes 
Asin U" = —c' (sin U—u’), since nuc=n'wc’. 
Alternatively this may be written: 

A sin U’=%' (sin U+sin I—sin I’—sin U’),  _......... (5) 
for, as a=1, uw’ =u+7—12'’=sin U+sin [—sin I’, where z and J represent the 
angles of incidence for the paraxial and marginal ray respectively, and 7’ and /’ 
the corresponding values for the refracted ray. 

The expression within the brackets will henceforth be designated as A(UU’). 
Its value for any given ray of semi-aperture U and U’ before and after refraction 
may be readily determined (by interpolation if necessary) from a set of specially 
prepared tables * in which the value of sin 7—sin J’ is tabulated for standard 
ee of the deviation D= U'’— U=I-—TI’, and for the relative refractive index | 
n' |n. 

Alternatively sin J—sin I’ may be replaced by i—7’=8, so that A(UU’) 
={(u+6)—sin (U+D)}, the value of D, in terms of 8 being given by the series 

Daven ate n'>—n?® hen? 
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For if we expand J in terms of sin J, and J’ in terms of sin I’, we have 
I=sin J—}.4 sin? 1+4.2.3 sin? J... . etc., 
and similarly for J’. 


* ‘The writer prepared some time ago a set of such tables, covering deviations between 0° and 12° 
and values of log n’/n between 18000 and 22000, 
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Hence 
, . 7” 3 _ 43 * G 5 
D=I-I' =(i-i )-1.4(S a) sete, since sin I’=" sin J and i=sin J. 
n 
. i _ 2 
Then equation (6) follows, since 6=i—7/ = Ailes i ae 
4 n’? n (n'—n)? 
= ———— 
(n’ —n)> 
The series in equation (6) is rapidly convergent for small values of 5, but for 
the rapid evaluation of A(UU’) it is quicker to pick out from a table of logarithmic 


sines a pair of values for log sin J and log sin J’, which differ by log n’/n, while 
I and I’ differ by the required deviation U’— U. 


adi 


6°; and so on. 


§4. AXIAL SPHERICAL ABERRATION OF THE MARGINAL RAY AFTER 
REFRACTION BY 4q-CENTRED SPHERICAL SURFACES 


We can now substitute in equation (4) the value of A found in equation (5) 
above, noting that the cy’ in the former equation corresponds to the c’ in equation 
(5) for each successive surface. 

Hence 


/ id 2 MY 4 pe 
Ngug sin U, A,=n,u, sin U,A,+ & ny'Uy'cy’ A(U,U,'). ..... (7) 
p=q 


In the absence of initial spherical aberration, the condition that the final 
ray U,,’ shall be free from spherical aberration is 


p=1 
fig (Ey Ge sN( Cla Cet) == Oat cu: a yng acs teeters (8) 
p=4q 


‘Thus the general procedure in building up a system free from axial spherical 
aberration is first to define the path of the marginal ray through the system by 
specifying the U’ of each successive surface (the U, of the pth surface being 
identical with the U;_, of the preceding). ‘These deviations having been chosen 
in accordance with the requirements of achromatism and the desired angular 
magnification of the system, factors are then selected representing the partial 
product 7,,'u,,'c,' for each surface, in such a way as to make the algebraic sum on 
the left-hand side of equation (7) zero. ‘This selection must of course be con- 
ditioned by other practical requirements of the system, such as the lens thick- 
nesses, the separations of the individual lenses, and the intersection heights, as 
well as by the need to remove other aberrations. 


§5. CONTROL OF COMA 


Control of coma will normally be achieved by the well-known application of 
the sine condition. With an initial u,=sin U, it will merely be necessary to see 
that final w,’=sin U,’. 

In the early stages of laying out the system, however, wu,’ has not been deter- 
mined, its exact value being dependent on the linear quantity c, and we only 
know the product wc. If, however, we write 
p21 aN COL CIB Ye 

sin U,’ 
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k= tunc,- and. Ky = 
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then the ratio of the magnification M, of the marginal ray to the magnification 


m, of the paraxial ray is 
hy(Ry +. Ky)(Ro+ Ko)(Rg t+ Ke). 2+ +--+ (ee (9) 
he(hy EER Ee ee (Rp+Kpa) 

which must be equal to unity if the final coma is to be zero. 


A 
(14 4,\(14 4s) 144s)... (14) 
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CAO) (ite) 
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(14 $2) (14 F) (14 32) sgl ee ee (1+) 

= 1 ky 3 ‘p 
pea. KG K ae 
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(14 aCe = i =) ( a) 


For 


which reduces to equation (9) when Ky and K, are each zero, which will always | 


be the case if the initial and final spherical aberration be zero. 


§6. SPECIAL CASE WHEN THE REFRACTING SURFACE IS PLANE 
OR OF EXCEPTIONALLY LARGE RADIUS OF CURVATURE 


When 7 is very large or infinite, c is also very large or infinite, while sin J is 
equal or nearly equal to —sin U and sin I’ to —sin U’. Hence A(UU’) tends to 


vanish identically. ‘The following substitute formula for nucA(UU’) may |} 


then be used :— 
Since nuc=nri, and sin J=c sin U/r, sin I’ =c’ sin U’/r, we have, substituting 
in equation (5), 
nuc(sin U+sin I —sin J’ —sin U’)=nri{sin U(1 +¢/r)—sin U'(1+¢'/r)} 
=n(Lsin U =D sin 0)" 9 eee (10) 
the rs cancelling out. 


This formula gives more accurate results for large radii and can be used for |} 


a plane surface. Otherwise it is less convenient as it involves two linear constants 
instead of one. 


To conclude this section it may be observed that the advantages in the use |} 


of equation (6) will be found more apparent in the case of separated compound 
lenses than for thin systems of conventional type, where the possibilities of 
combination of glasses and surfaces have been thoroughly explored. 


Even where trigonometrical methods are preferred, however, the quantity | 


nucA(UU") should still be formed for each surface parallel to the other work, 


since it is so readily calculated from the available elements. It shows at once |} 


how much each surface is contributing to the total spherical aberration, reveals 
the necessity for adjustments at an earlier stage, and shows to what extent the 


spherical aberration has accumulated beyond the power of subsequent surfaces 
to reverse. 


§7. CURVATURE OF FIELD FOR CENTRAL RAYS 


The formation of images by thin pencils passing through the centre of curva- | 


ture of each refracting surface can seldom be realized in practice. "The Petzval 
surface on which such theoretical images would be formed is, however, of con- 
siderable practical importance, since the distance of the real intersections of the 
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image-forming rays, before or behind this surface, is the measure of the degree 
of astigmatism introduced by the refracting system, and, as is well known, the 
formation of a plane anastigmatic image of a plane object requires as a first 
condition that the Petzval curvature shall also be zero. 

‘The Petzval curvature for indefinitely thin pencils of central rays may be 
defined by the well-known formula 


f Daly! = 
Zia _ 2x, 45 Mp — 
/ 12 2 / , 
iN Bikes genta re 


where 7, is the radius of curvature of the pth refracting surface, x, and Y, are the 
rectangular coordinates of an extra-axial object point, referred to an origin B 
on the optical axis where it is cut by the object surface, and «,' and Y,,’ are the 
corresponding rectangular coordinates of the pth image of the point x, Y, referred 
to the pth paraxial image of B as origin. 

This equation is an approximation only, and although for most conjugate 
positions of object and image and a fairly small angle of field it gives sufficiently 
accurate results, with the development of wide-angle systems and of less con- 
ventional types there would appear to be good grounds for replacing, wherever 
possible, the approximate Seidel conditions with exact formulae. Coma is 
already taken care of by the sine condition; an exact addition equation for axial 
spherical aberration has been supplied in a previous section; summations are 
now given for the Petzval and sagittal curvatures. 

An exact form of the Petzval equation for a single surface follows from equation 


(25), viz., 
Rie aA (3) 


For, following the well known method for finding the Petzval equation, with 
B as origin, let x, Y (figure 2) be the coordinates of an extra-axial point in the 


(e) hi xth oh x-h 


—xN'5 e— r*—— 


Figure 2. 


object surface and x’, Y’ the coordinates of the corresponding extra-axial image 
point formed by an indefinitely thin pencil of central rays. With centre at the 
centre of curvature of the refracting surface, describe arcs through the points 
‘xY and x'Y’. Let h and h’ be the respective depths of curvature of these two 
arcs. We then have 

Nae ho A= x — hy | Vom. 
Substituting in the above equation for A/’ and Al, we have 


pag at peo Ge **). 
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Dividing by n' Y’=n'm?Y? and transposing, 
x x >( 1 1 )-m (x—h)? 
&. oe ngellmrape © Saye 0 


WY? AY? OY \wv ey os es 
- hlY=h'/Y’. 
since / p V24h 2Ch 
Remembering that Y’/Y=C’/C and C"/Y = C/Y, and that —— aye cyan 
x’ Cane elite he _ot («— i). 
n VY? nY2 2 nn'r (1 7 y2 OE: 


Again, if 8 is the semi-angle of field, iss tan 44, so that, finally, 


iS EL! ln’ — ma oe n’ ce a 
pv ave Oe 1+tan? 0-2  mo($ tan 10) }. = (11) 


Since the «’/n'Y” of one surface corresponds to the «/nY? of the next, the 


exact position of any image point formed by indefinitely thin central pencils in 
the Petzval surface after g refractions is given by the formula 


ay = p= i 12-9 a 
Yee | se {1 +-tan?} 56, -22 "p hy (3 —tan $0 lee 
Mp 
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If the initial object surface is plane then is Zero. 


er 
mY, 
‘The right-hand side of aa (11) then becomes 


es 24 n' 


(a) For a plane refracting surface, 1.e. r= ©. 


‘This can vanish: 


(>) When the object is at the centre of curvature of the refracting surface, 
i.e. when m=n/n’ and tand0=1. _ 

Combining these two conditions, it will be seen that the virtual image of a 
plane object formed by a plano-convex lens of focal length F at a distance 
(n' —n)F/n' from the object has no Petzval curvature, although by the approximate 
formula the curvature is 1/nF. ‘This relative position of object and image is 
approximately realized at the field lens of a Huygenian eyepiece and in the front 
combination of some microscope objectives, thus assisting in the flattening of 
the field in such cases. Unfortunately, in a cemented plano-convex achromatic 
lens, the tendency of the crown is to neutralize this advantage, unless the refrac- 
tive index of the glass of lower dispersion is at least as high as that of the flint. 

Equation (11) consists of a principal term representing the usual Petzval 
sum with the addition of a corrective term for each surface. However, the exact 


condition that the Petzval sum may be zero can be expressed in a simpler form ~ 


by dividing the axial magnification formula (1) by Y’ instead of by Y” and 
replacing Al and Al’ by (~—h) and («’ —h’) as before. 
Then for a single surface we have 


x x 
Ue sg ae us = us — Stan S050 - ot Eee (13) 
where 6 is the paraxial deviation for a ray of semi-aperture u and wu’ before and 


after refraction. 
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For a succession of spherical surfaces the ane becomes 


Ua x p= 
no = mg — ES Op tds Op, sas nw bvists ak (14) 
so that if the object and final i image surfaces are each to be plane, 
p=1 
200, cat S0 == OM 0 Te seteve (14a) 


p=4 
§8. ADDITION FORMULA FOR SAGITTAL ASTIGMATISM 

By using equation (3) a condition similar to equation (13) above may be 
obtained for the curvature of the sagittal field. 

Dividing by n’Y’ sin V’=nY sin V, and again substituting (X’—H’) for 
AL’ and (X —H) for AL, and transposing, we have 

es = pes (i: Se se ARCENASS 
Na NA we nY sin V 
n'u'c' A(VV’) 
=a ee (15) 
using capital X and capital H, V, V’ and Vin place of x, h, u, u’ and @ to distinguish 
the sagittal ray. 
Then for a succession of g surfaces 
2G DC a Bs SIN GA ED 
Ti Gps ag = OE eo See A Siesta 

The final term in equation (16) has one serious disadvantage. It has already 
been stated that in the spherical-aberration formula (5), the scale of each element 
of the complete system is determined by a single linear constant c’. ‘This.constant 
again appears in equation (16), but there is no longer complete freedom of 
selection. 

Although the scale of the system as a whole is still open, the relative scale as 
between the separate components has already been pre-determined by the 
selection (in conjunction with the paraxial deviation 5) of the angular subtense 
WY, of the field presented to each surface, thus fixing the positions of the centres 
of curvature and leaving no remaining latitude in the radii after these have 
produced the prescribed deviations. 

It is therefore desirable to expand this term as follows :— 

n'u'c A(VV") _ nucA(VV") _ wav) (4) 
We sian wy sin V sin V i 
a rue (cosec OY 4 = —tan ry). ipa (17) 


sin 


=u — Stan} + 


A(VV’") cosec 


For the preliminary work it will be sufficient to take nV’ 


as the approximate value of the above term. The remaining corrections can be 
applied at a later stage. 

Any combination of the second and third terms on the right-hand side of 
equation (16) will give a flat field for sagittal rays provided the total sum is zero; 
but to ensure that the sagittal image is formed in the Petzval surface, equation 
(14) must also be satisfied; or, a little less precisely, the two sums in equation 
(16) must be separately zero. 
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MODULUS FOR SHORT WIRES 1) 
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ABSTRACT. The method consists essentially in the measurement of the time for al 
quarter vibration of a short wire clamped at one end. The wire breaks a circuit as it 
starts its transverse vibration, and is made to. strike a very light lever and open another 
circuit just as it completes one-quarter of its first vibration. ‘The time between the two 
events is measured ballistically. From a knowledge of the time of vibration, the length, 
the diameter and the density of the wire, Young’s modulus can be readily calculated }} 
The method is particularly suitable for thin wires having a diameter of the order of 
0-015 cm. and a length of about 2cm. The results are consistent to a few per cent. 
Comparison of the results with the accepted values for five different wires shows a good 
agreement. 


$1 IN TRODUCALTON 
lo: RAYLEIGH (1926) correlated the frequency of transverse vibrations 


of a bar with its Young’s modulus, and the correlation was used, pro-f 

bably for the first time, by Prosad (1929), who maintained the vibrations 
electrically, and determined the frequency chronographically by connecting to 
the free end an attachment for which an end-correction was applied. ‘This 
correction amounted to more than 200% for lengths of the order of several cm. 
‘he mean results of Prosad differ by about 5% from the results he obtained b 
the flexural method, while his curves suggest that a divergence of about 15%] 
from the average is possible. 


In the present work, the frequency is calculated from the time of one-quarte 
of a vibration of the wire. 


§2. EXPERIMENTAL DETAILS 


A length of a few cm. of a thin straight wire is used. One end of the wire} 
is rigidly clamped, and the other end can execute lateral vibrations. It is} 
arranged that as the free end of the wire is released from a depressed position\ 
it breaks a circuit, and just as it reaches its normal position it opens another} 
circuit. ‘The time interval between the two events is measured by a condense | 
circuit described by Klopsteg (1920). : 

The principle will be understood by reference to figure 1. The two contacts} 
K, and Ky are opened in quick succession at the beginning and the end of thel} 
short interval to be measured. Then the battery is reversed, K, short-circuited] 
by Kg, and the two-way key K; put to the galvanometer side. If the setting off 


the resistances R and S has been right, the galvanometer throw will be zero, and} 
the interval will be given by 


t=CRlog, R/S. 
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The contact K,. ‘The wire X under test is clamped at C (figure 2) by means of 
two plane-edged aluminium jaws which can be tightened on the wire by nut and 
screw. ‘T’he rigid clamp forms one terminal A of the contact. ‘The clamped 
end of the wire is sharply defined by the plane edge of the clamp. ‘The free end 


Figure 1. The principle of the method. 


K of the wire can be depressed and kept in a constrained position by means of a 
stout nickel wire D having a small, flat, horizontal edge, with which it makes 
good electrical contact. D is connected to the second terminal B. 

The contact at K is made at the extreme edges of X and D, so that when D 
is suddenly pulled outwards by means of an electromagnet or an attached thread, 


Figure 2. he contact Ky. 


the break of the battery-circuit of figure 1 is abrupt and the current drops 
immediately to zero. 

The position of D can be adjusted both horizontally and vertically by means 
of the screws E and F, and thus measurements can be carried out for various 


lengths of X and for various amplitudes of vibration. 
The arrangement is fitted on an ebonite base G, which itself is fixed centrally 


on a levelling table. 

The contact K,. The contact K, has a very small inertia and at the same 
time has a reasonably low and fairly constant resistance (0-7 to 0-9 ohms). 

JOL (figure 3) is a fine steel wire (diam. =0-02 cm., length =5 cm.) pivoted 


414 Y, Ln ousep 
at O, where a pin passes through a small loop in the wire. The pin is soldered tay | 
a copper plate P connected to a terminal screw Ty. The end L is bent at righ | | 
angles to the length JL, and a short thin platinum wire (diam. =0-0075 cm. }} ? 
is soldered at N, where the lever rests on the plane top of the copper wire Sj] 

which is freshly but thinly amalgamated with mercury, and which forms the second | 


Figure 3. The contact Ky. 


lead. A very thin, loose wire W of copper (S.W.G. 40) is soldered between 
JL and the plate P to improve the contact between the lever and the terminal T,-. 
It also helps, by the slight pressure it exerts, in keeping N and S in good contact, 
and thus the resistance between T, and T, is small. ‘The whole arrangement is 
mounted on an ebonite base E, which is held by a fixed stand. 


Figure 4. The experimental arrangement. 


The vibrating wire X (figure 2) passes between E and JL, and is arranged} 
to strike the lever normally at some point Q near the end L. As the lever receives] 
the impulse, the contact between N and S is suddenly broken, and the point L} 
is raised enough to come into the field of action of the small magnet M, to which }] 
it is attracted, and thus the contact remains open. | 

A force of a few centigrammes weight is sufficient to pull off the point N] 
from S. 

The actual arrangement. For Ks, Ky and K, of figure 1 we used mercury 


keys of the rocking type. ‘They are constructed on one block of parafiin wax, , 
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nd are all operated by a common handle H (figure 4). The lengths of the 
10vable links between the mercury pools are such that on switching over from 
he first position to the second, the short-circuiting key K, and the reversing 
ey K, are closed just before the discharge K,. 
R and S are two large resistance boxes ; C is a standard mica condenser of 
apacity 2:0005 x 10-* F. ; G is a ballistic galvanometer having a sensitivity of 
scale divisions per micro-coulomb, and a deflection of about 1/10 of a scale 
ivision can be detected ; the battery B has a constant e.m.f. of about 32 volts. 


Adjustment of the test wire. ‘The position of the unconstrained test wire, 
rhich is situated normally with its free end below the end of the lever, is carefully 
djusted by levelling until it starts to raise the lever, and the contact K, is just 
roken. ‘The break may be judged by viewing the tip N and its image on the 
right amalgamated surface of S through a ens (figure 3). A white background 
ssists in producing a distinct view. ‘The right position may be located to within 
nuch less than corresponds to a 20° rotation by one of the screws of the levelling 
able supporting K,. ‘This rotation is equivalent to a linear displacement of less 
han 1/200 of the smallest amplitude of vibration we used, which is about 3 mm, 

On account of the very small inertia of the lever, the test wire is not appre- 
iably bent as it starts to raise the lever. ‘This is especially so for short lengths 
f the test wire, since the bending varies directly as the cube of the length. 

When the “zero” position of the test wire is exactly located, it is depressed 
o the point D (figure 2) to close the contact K,. Care is always taken that the 
lepression of the free end does not exceed 1/6 of the length used. For this depres- 
ion, the initial amplitude of vibration may be regarded as small, and the wire will 
ot change its zero position. 

By virtue of the large velocity of vibration near the mean position, an error of 
°/, of the amplitude arising from a slight error in locating the test wire will pro- 


2 ; ree a 
luce one of only—z % in the time of one-quarter vibration, or an error of 8 % 
7 
n Young’s modulus. 


§3. THE METHOD OF MEASUREMENT 

With K, closed (figure 4), observations are made by first approximately 
djusting the position of the test wire, and then closing K, and taking trial values 
or Rand S. On releasing the test wire, the contacts K, and K, will be opened 
1 succession after a quarter vibration. ‘The keys K3, K, and K, respectively 
re then immediately switched over by the handle H. If G shows a deflection, 
eadjustment of the resistances R and S is made until the deflection is reduced 
ractically to zero. It is generally sufficient to change one of the resistances only, 
nd keep the other constant at a suitable value. ‘The test wire is then exactly 
ycated, and the final adjustment of the resistance is made to give zero deflection 
xactly. 

The values of R and S are now recorded, and a new balance is obtained with 
ther values of the resistances, the zero position of the test wire being always 
hecked during the final adjustments. An average time can thus be obtained. 

If the contact K, has a small resistance 7, then 


t = C(R+r)log, a 
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The frequency of transverse vibrations will be given by 
Ei 
Se | 
Lord Rayleigh (loc. cit.) has shown that, for lateral vibrations of a bar 0) 
length /, density p and Young’s eee Y, the frequency is given by | 
km 
ar, /(Y/p), 
where k is the radius of gyration of the section about an axis perpendicular to thd | 
plane of bending, and m is an abstract number whose value for the gravest mods 
of vibration of a clamped-free rod is 1-8751. | 
For a circular wire of radius a, k= 4a, and hence 


ye Tea | 
s | 


from which it will be seen that for a given wire v*/* is a constant quantity. 

The length of the wire between the clamped and the free end is measured by 
callipers, and its average diameter is measured by a micrometer screw gauge.|} 
The density is determined from the mass and the dimensions of the sample | 
frequency determinations are made for various lengths of the wire, and the mear) 
value of v?/4 is obtained and used for calculating Y. 


§4. RESULTS 


Typical results are given in tables 1, 2 and 3, where the c.g.s. units areé 
used throughout. The data in table 1 show the accuracy with which thd 


Table 1. Accuracy in measuring the time 


R (ohms) S (ohms) . t (sec.) Average t % deviation 
1500 699 0:00229 —1:-4 
1250 514 f 0-:00222 alates 

0-00226 
1000 322 0-00227 —0-4 


750 168 0:00225 +0-4 | 


Table 2. Effect of change of amplitude 


All R (ohms) S (ohms) t (sec.) Mean t 
2000 1010 0:002734 
0-222 0-:002738 
1500 602 0-002741 | 
2000 1013 0-002723 
0-109 0-002715 
1500 609 0-002707 
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adjustments can be made and the time measured. The observations were made 
on a constantan wire of length 5-535 cm. and radius 0-0273 cm. 

The maximum departure from the average is less than 2%. 

The figures in table 2 show the effect of change of amplitude on the time off : 
quarter vibration. ‘The observations were made on a nichrome wire of length} 
5-13 cm. and radius 0:01812 cm. The times for two widely different ratios of} 
amplitude (A) to length (/) are given. 

It will be seen that for the ratios of A// given, the effect of change of ampliendd 
does not exceed 1%, which is of the same order of magnitude as the experimental }f 
error. However, the ratios of A/I we actually used in our determinations were 
always less than 0-16. | 

A summary of typical results at room temperature (18° c.) for six different! 
wires is given in table 3, where Y, stands for Young’s modulus by the author’ si 
method, and Y, stands for the accepted value (Kaye and Laby, 1928; Childs, J 
1943). ) 

The table shows how far the results of different observations on the same} 
wire are concordant among themselves. It will be also noticed that, in general | 
there is a good agreement between the observed and the accepted values of the) 
elasticity. 


§5. CONCLUSION 


When measures are taken to secure good clamping of the test wire by a massive } 
support, small inertia of the lever, high electrical conductivity at the contacts 
K, and K,, small area of contact with the tip of the lever, exact location of the test 
wire, good insulation of the various keys, and sudden break of the battery circuit } 
by K,, the method appears to be very useful in measuring, with ordinary laboratory 
equipment, Young’s modulus of elasticity, particularly for thin and short wires. 
No correction is required, and the wire need not be especially shaped or adapted. 
The average result will give a good measure of the true modulus of elasticity. 


I wish to express my thanks to Dr. M. A. El-Sherbini for his valuable help: 
and advice during the progress of this work. 
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ABSTRACT. A method is described which permits the computation of the contributions 
made by the individual surfaces of a lens system to the tangential aberrations of the final 
image for a pencil of finite aperture and obliquity. It employs the transfer coefficients 
previously described and forms part of a general development of the trigonometrical 
method of design. When applied to pencils of small aperture and small obliquity it gives 
an analysis of the surface contributions in exact agreement with known primary methods, 
and thus provides an extension from the analysis of the primary to the trigonometrically 
determined aberrations. In particular, expressions are derived for the contributions to 
the chromatic aberrations made by any component of the system which are simple functions 
of the dispersion of the component. 


Si IN TRODUCTLON 


value to know the contributions which each surface makes to the image 

aberrations. Methods of computing the surface contributions to the 
primary aberrations are well known, but, for aberrations other than these, no 
eneral analysis into surface contributions is available. Using the properties of 
the differential transfer coefficients developed in Part I of this paper, it 1s possible 
0 solve this problem for the tangential aberrations of a system. We proceed to 
develop a trigonometrical method of computing the surface contributions char- 
icteristic of a lens system for a given aperture and position in the image field. 


ff the course of the design of an optical system it is frequently of considerable 


§2. THE CONTRIBUTIONS TO THE LONGITUDINAL AND 
TRANSVERSE CHROMATIC ABERRATIONS 


In Part I an expression is derived for the longitudinal chromatic aberration 
»f the system at a given zone in the form 


aL! 

(pas aa ve St a ee 67 
Efe adie Ne Te To IN 68 

pcre aN, Om Nyy Om N,2’ CG 


he summation being extended over all the components of the system. Each 
erm in the summation in (67) represents the contribution made by a component 
o the longitudinal chromatic aberration of the final image. ‘Thus for the general 
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component, h, the trigonometrically determined contribution to the longitudinal] 
chromatic aberration, which we denote by T'LchC,,’, is given by 


oa (NON), ee (69) 
If the surface contributions are required rather than the component contributions,} 
it is obvious that they may be easily obtained. Each term in the summation off} 
equation (69) is derived from a pair of terms as seen in equation (68), the twa 
members of each pair referring to single surfaces. Hence the longitudinal chro- 
matic aberration may be expressed as a summation over the individual surfaces 
of the system if desired. For each air-glass surface there is a single term, while 
for a cemented surface there are two-terms each of which involves the partial 
dispersion of one of the glasses forming the contact surface. ‘The term, or the 
sum of the two terms, relating to each surface gives the contribution of that 
surface to the final value of Lch’. For the paraxial chromatic aberration we 
have similarly 


T LHe = 


: 5 a’ 
WBE a in (N,—N,),,  ° > = eee (70) 
; = 
whence TIlchC ens dV ne Se (71) 


In Part I an expression was also Pee for the transverse chromatic aberration} 
for a given obliquity in the form 


amas OH bra ms 3 
Tch' => | ae ie Aas. hee (72) 


This expression is based on the trace of a principal ray of the desired obliquity} 
using the refractive indices, of the intermediate colour d. The subscript L 
indicates that the plane in which the ordinates, H’, are measured is a fixed planeJ}} 
It is convenient to use the paraxial image plane for this purpose. In a manne | 
similar to that followed in the case of the longitudinal aberration, it may be show ! 
that the trigonometrically determined contribution to the transverse chromatic 


aberration for the general component, h, is given by 


TTchC,! = | S| (NAN eo es (73) 
h LT’ 


If the surface contributions are desired, the equation (72) may be rearranged as} 
a summation over the surfaces in a manner similar to that already discussed in} 
the case of the longitudinal aberration- 

The contributions made by the components of the system to the final]| 
chromatic aberrations are frequently of great value in practical designing, being} 
much more useful than the contributions made by the surface of the system, 
Their value lies in the fact that they are such simple functions of the dispersions, |} 
and it is often possible to improve the achromatism of the system by a change of] 
glass at one or more components after inspection of the values of the contri-} 
butions. 

In the accompanying table, the results of the computation of the contributions] 
made by the surfaces of a photographic objective to the paraxial longitudinal} 


Transfer coefficients for use in the design of lens systems: III 421 


chromatic aberration and to the transverse chromatic aberration at two obliquities 
are compared with the corresponding values for the contributions to the primary 
chromatic aberrations as calculated by the method of Conrady (1929). 


Table 
Conrady Conrady Conrady 
Surface TlchC’ IchC’ ch Cx TchC’ EL ches Lche- 
Bh 10° 

1 3:5406 3-5406 —0-0086 —0-0083 —(0-0223 —0-0210 
D 3°8759 3-8759 0:0778 0:0767 0:2138 0:1931 
3 —6:6646 —6-6646 —0-1331 —0-1313 —0-3659 —():3304 
+ 0:4251 0-4251 0:0723 0-0713 0-1949 0-1795 
5 —2-0300 —2-0300 0-0101 0-0099 0:0271 0-0250 
6 0-0052 0-0052 —0-0124 —0-0122 —0-0341 —0-0308 
7: 1:6591 1:6591 0:0034 0:0034 0:0089 0-0086 


In the first place it will be seen that the contributions determined by the two 
methods lead to identical results for the longitudinal aberration. Secondly, 
the table shows that at an obliquity of 4° the values given by the two methods are 
in fairly close agreement, the small difference being due to aberrations of higher 
order than the primary, which are already present. As the obliquity increases, 
the two sets of values diverge further, as is to be expected. At very small obliquities 
the two sets of values would be in exact agreement. It will thus be seen that a 
perfectly general trigonometrical method of calculating the surface and component 
contributions for the chromatic aberrations of a system has been established. It 
is accurate, as may be checked by full ray-traces, and when applied to the paraxial 
region and to pencils of small obliquity it gives results in exact agreement with 
existing primary theory, thereby becoming a satisfactory extension of known 
methods. 


§3. THE SURFACE CONTRIBUTIONS TO THE SPHERICAL ABERRATION 


In the ray-trace of the axial pencil we have traced rays incident at certain 
zones of the lens system. In the image formed after refraction at surface 7 of 
the system, there is, for any such zone, a spherical aberration given by 


LA; =1/ — Ly’, 
which may be expressed in angular measure as 
AA, =LA; sin U;'/S, 
eioA sin U cos Us /(iy —X,).  nexesn (74) 


The angle AA, measures the departure of the ray from the ideal direction after 
refraction at the surface z. If the ray could be turned through an angle dU;' = 
—AAj’, the spherical aberration behind the surface z would be reduced to zero, 
ind the effect at the final image would be the introduction of an amount of spherical 
iberration opposite to that introduced by the actual refractions over the first 7 
surfaces of the system. In other words, the shift of the intersection point of 


he axial ray with the principal axis in the final image space due to this imagined 
28-2 


422 F.. D. Cruickshank 


rotation of the ray through the angle dU,’ = — AA;’, at surface 7, will provide a | 
measure of the sum of the spherical contributions of the first z surfaces of the |}: 


system. We shall denote this quantity by the symbol & TSC’, the letters standing | | 
for the trigonometrically determined spherical contribution. Hence | 


¢ OL. / OL. , 
‘ = — AAS op Mate rf 
bap pele = BU eS aU; AA;, (75) 
and the individual contributions from the es are given by 
TSC ga y TSC’ — Ss TSG.) ‘> Soe eee (76) 


This gives a general trigonometrical method of analysing the surface contribu- 
tions to the final spherical aberration. Little additional calculation is involved | 
when the general methods outlined in these papers are employed, for the necessary | 
transfer coefficients are already obtained in the main computation. ‘The only 
limitation to the accuracy of the method is that differential transfer coefficients 
are used to calculate the effects in the final image-space of the rotations, dU’, | 
of the ray, If the latter are large, the transfer coefficients will not predict their. 
effects accurately. This difficulty is met to a very large extent by using appro- 
priate second-order correction terms, as developed in the previous paper. If ] 
dU,,’ is the change in the direction of the final emergent-ray due to the rotation } 
of the ray through dU,’ behind surface 7, then the shift of the intersection point 
along the principal axis is by equation (46) : 


él,’ = C(U,’) cosec U;, dU;(1—cot U,,/ dU,,'), ai ovtexieneis (77) 
and with due attention to the sign of the longitudinal aberration this leads to 
TSC =C(U;) cosec VU, dU{(V—cot U, dU een (78) 


This equation (78), incorporating the correction term, should take precedence 


over the earlier equation (75) and be used for the normal calculation of : TSC% 
The separate surface contributions then follow from equation (76). 


§4. THE SUREACE CONTRIBUTIONS LO THE DISTORTION 


From the results of the trace of an oblique pencil through the lens system we } 
can obtain the value of the distortion for this obliquity in the images formed after 
refraction at the successive surfaces of the system. In the usual linear measure 


this is 
(dist’); =(Hia— Hyrs)is 
and in angular measure 


AD,’ = dist,’ cos U4 Syrsi 

= dist, cos* U7, nj (LX 5). eee eee (79) 
This angle, AD,’ measures the departure of the principal ray from the ideal-|} 
direction, as regards distortion, after refraction at surface 7. If we could swing | 
the principal ray through an angle dU;,,;= —AD,’, we should restore the ray to | 
its ideal direction and thereby remove from the final image the amount of distortion | 
introduced by the refractions at the first 7 surfaces of the system. The effect 
of such an imagined rotation then is to introduce into the final image an amount | 
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of distortion opposite to that contributed by the first 7 surfaces of the system, 


an amount which we denote by = TDC’. Hence we have 


-UTDC'=. d(Hi,—Hi,,,)= —dH 


prf 
= 0 OF rs ’ 0 ee , 
oes Te LU i = a dU, AD, 
. y 0H, 
that is, 2 TUG = =P Al)! 8 
SS amare (80) 


t 
and the contribution from the surface z is ues by 


TDGSSTDC- STDC (81) 


Equation (80) is adequate provided the angles dU/,,; are small. _ Improved 
accuracy is obtained by using the second-order correction term ; the difference 
in the amount of computation involved is so little that the corrected equation 
(82) given below should be used as a matter of course. Using equation (49) 
of the previous paper, and omitting the subscript pr, as all the symbols refer 
to the principal ray, we obtain finally 


UD Cl Clr) psecUe alii (lta Ur dg)e =. ae (82) 


§5. THE SURFACE CONTRIBUTIONS TO THE TANGENTIAL COMA 


We begin by computing the value of the tangential coma for a pencil of given 
obliquity in the images formed after refraction at the successive surfaces of the 
system. If we choose to regard the a and 6 rays of the pencil as aberrant, then 
it is the displacement of their intersection point, Ab, from the intersection point 
Pr which afflicts the system with coma. If after refraction at surface 7 we could 
swing the a and 6 rays until both passed through the point Pr, the tangential 
coma originally present in the image behind surface 7 would be removed. Such 
a swing of the rays would also leave the curvature of the tangential field unchanged. 
The effect of such an imagined rotation of the rays on the final image would be 
the introduction of an amount of tangential coma opposite to that contributed 
by the first 2 surfaces of the system, that is, an amount which we denote by 


— TCC’. Thus we have a basis for the derivation of equations for computing 
the contribution of each surface of the system to the ee coma of the 
final image. 

The angle through which the a ray must be rotated is 


a iy / , 
dUj,;= —comar, cos Uj,/'S); 


= —comay,cos* Uj,/(Lip—Xq )ie tt ee ee (83) 
The corresponding angle for the rotation of the 6 ray is 
dU;, = —comar, cos? U;/(Li,—Xp i we eee (84) 


In order to calculate the contributions at the final image we must follow the shifts 
of the intersection pointsAb and Pr in the final image space due to these rotations. 


The shifts are 


Bie ena eet 


eearg iii if «| a enere 85 
aU%, erg ss! eS 
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OMe guy, + a ay; 86 

dH iy = ay Ou dU ye au}, bi» c ( ) | | 

dH, = dL, tan. Cork 0) eae (87) | 

It then follows that ; 
SCC =dH,,,-diy,, ~ 2 eee (88) 


and the contribution made by the surface 7 1 is ae by | 


TCC) =S37CC 2 370C) (89) 


As in the previous case, it is generally advisable to use the second-order correction |}} 


terms. It will be easily seen that the modifications necessary in equations 
(85) to (87) are that the first term on the right-hand side of (85) and (86) is multi- | 
plied by [1 —cot (U,'—U,’),d@U;,], and the second term by the corresponding 
quantity [1—cot(U,’ — U,’),,dU4,]. 

The essence of the foregoing method is the computation of the effect at the) 
final image of a change made in the refracted rays behind surface z, the change 
being such that the coma present behind this surface is removed. Another) 
change which would remove the coma in the image behind the surface z is a swing | 


of the principal ray until it passes through the point Ab,. It may seem more} 
logical to make the changes in the directions of the rays a and b, because we regard |} 


coma as due to the aberrant behaviour of the outer rays of the pencil, but it is 


simpler to make the change of the principal ray for the purpose of calculating |} 
the coma contributions, as it involves less than half the amount of computation. jf} 


The fact that the contributions computed either way are in close agreement 


confirms the validity of the general method. ‘The angle through which the }ff 


principal ray must be rotated to remove the coma behind surface 7 is 


Jf £ 
dU), ,,= coma; cos U;, fae pri 


=coma7; cos? U7 /(L,— Xone  eee (90) 
Following out the effect in the final image space of such a rotation we have 
: Re OL ees 
CC = dU, Vor is (91) 


An additional advantage of the second method is the ease with which the cor- 


rection term may be applied. As in the case of distortion, we obtain 
TCC = — CU ;)p, sec Uppy dU: +tan U;'dU,,,). ...2.- (92) 


Finally, then, we use equations (90), (92) and (89) for the computation of the} 
contributions to the tangential coma. 


§6. THE SURFACE CONTRIBUTIONS TO THE CURVATURE OF 
THE TANGENTIAL FIELD 


The curvature of the tangential field is measured usually by the distance of ]} 


the intersection point Ab from a plane at right angles to the axis through the 


focus of some chosen axial ray. his ray is usually the extreme ray of an axial |} 


pencil having the same relative aperture as. the oblique pencil. ‘Thus 


: | Xq’ = Lay — Ln» |i 


| 
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As long as the point Ab lies in the plane through the intersection point of ‘the 
ray m with the principal axis, the curvature of the tangential field would be zero 
according to the last equation, but the coma would be profoundly affected by 
the H’ co-ordinate of the point. Considering the curvature as quite distinct 
from the coma, we could say that in a comatic system the ideal location of the 
point Ab as far as curvature is concerned is in the plane of the m focus at a point 
distant (coma), ‘from the intersection point; Prm, of the principal ray with this 
plane. ‘This analysis of coma and curvature may appear somewhat arbitrary, 
but so are most measures and analyses of the aberrations of asystem. We propose 
to use the above ideas as a-basis for determining the surface contributions to the 
tangential curvature. 

From the results of the ray trace we calculate the tangential curvature of the 
image formed after refraction at successive surfaces of the system. Considering 
the rays a and b after refraction at surface 7, a rotation of each ray so that their 
intersection point Ab, now lies in the plane of the m-focus at a point distant 
(coma)7;; from the point Prm; would remove the curvature originally present 
in the image behind the surface 7. The effect of such imagined rotations would 
be to introduce at the final image an amount of curvature opposite to that intro- 


v 
duced by the first z surfaces of the system, that is, an amount —XU7TXC". ‘This 
leads, as before, to a measure of the surface contributions. 
The angle through which the a ray must be rotated is 


dU, =X 7(tan Uj;—tan Uj,,;) cos Uj,/\Si,; 


=A 7 A(tan Ue tary Ur.) cos sO (L a aNg pe esse (93) 
with a corresponding expression for the rotation of the b ray given by 
dU,,=X7(tan U;,,;—tan Uj,) cos? Uy (Lm —X5');- oe (94) 


The effect of such changes in the directions of the rays on the curvature of the 
final image is given by 


NoMa) ANB 
ea (Sati toes ah ge 95 
aX _' = dLiy = OU’. oe au, + aU! Ue bi ( ) 
and then SNC eet ey ns 7 eee (96) 
D i—1 
and ERC LON TN COS DX (eta ee dee (97) 


The correction term to be applied to equation (95) is exactly similar to that for 
equation (85), already noted. 
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ABSTRACT. A knowledge of the transfer coefficients for the aberrations permits the 
changes inthe residual aberrations of the system to be calculated for any small departure | 
from specified dimensions. Provided the changes in the residuals are small in com= 
parison with the residuals themselves, the effect of the departure from the specified 
dimensions will be negligible. Hence the calculated values of the transfer coefficients 
provide a basis for the estimation of the limits within which each dimension of an optical 
system must be controlled in the workroom if systems having uniform characteristics are 
to be produced. 


§1. INTRODUCTION 


duced by failure to achieve the curvatures and other dimensions specified 

in the design. On the other hand, the cost of production may be increased 
considerably by striving to attain the specified dimensions within limits which 
are unnecessarily fine. It is thus an important part of the designer’s work 
to specify the limits within which each quantity must be controlled during 
production. It is proposed now to examine this problem in the light of the 
analysis developed in the preceding papers of this series, and it will be shown 
that a method of estimating tolerances of this kind may be obtained which links 
on logically and conveniently with the general method of final design which 
has been developed. _ It is to be emphasized that the tolerances under considera- 
tion are not those which are generally termed ‘optical tolerances”’, which have 
relation to the amounts of residual aberrations which may be permitted in a-] 
system which aims at a certain standard of definition in the image. The present 
considerations relate to the degree of control to be exercised in the optical shop 
in the production of a given system if a uniform product is to be produced. 

_ The final design of a system always represents more or less of a compromise. — 
It is characterized by certain residual aberrations which have been calculated | 
during the final stages of the design, and which experience, or the performance }} 
of a carefully built prototype, shows to be compatible with satisfactory image ]| 
definition. ‘The whole success of known methods of design depends on the fact 
that the mathematical analysis of the aberrations gives a fairly reliable guide to 
the actual physical aberrations of the system when constructed. In the course 


I’ the production of an optical system, serious imperfections may be intro- 


Transfer coefficients for use in the design of lens systems: IV 427 


of production, any variation from the specified dimensions will result in a variation 
n the actual residual aberrations, but production will be satisfactory provided 
these variations are small compared with the residual aberrations themselves. 
Hence it is reasonable to base a system of estimating tolerances on the calculation 
of the effect on the residuals of a departure from the specified value of each 
quantity in the system. ‘he reliability of these tolerances will be of the same 
order as the reliability of the design methods as a whole. The system of transfer 
coefficients developed in Part I provides the machinery for the method, little 
additional computation being required. 


§2. THE TOLERANCES FOR THE CURVATURE OF 
EACH SURFACE OF THE SYSTEM 

In Part I it is shown that it is possible to calculate a transfer coefficient which 
measures the rate of change of any aberration with the curvature of any surface 
of the system. ‘These coefficients are calculated normally for the purpose of 
the final differential correction of the system. ‘They are now available for the 
further purpose of estimating the tolerances permissible during production. 

If we use the symbol A’ to denote any of the tangential aberrations of the 
system, then the calculated values of 0A’/dc reveal the effect of small changes of 
curvature on this aberration. After a close scrutiny of these quantities, tolerances 
may be set by assigning a permissible curvature variation at each surface such 
that the sum over all the surfaces of the effects on the residual aberrations due to 
the occurrance of such variations shall not exceeed some specified fractions of the 
values of the residual aberrations. If a curvature tolerance, +dc;, is selected 
at each surface in this way, the corresponding tolerance in the radius of curvature 
is +7;2dc;. It is often convenient to express the tolerance in terms of the number 
of fringes across a surface of a certain diameter. It is easily shown that the 
difference in curvature, dc, between two spherical surfaces of curvatures c and 
¢ + dc, which show «x fringes across a surface of diameter 2a when placed together, 
is given very closely by 


r 
de = atl Sse mare Me | nidcags (98) 
Writing 0A’ /dc, for the change in the aberration A’ per fringe change in curvature, 
we have for any surface 


“= 6 Seve) pe a es (99) 


Equation (99) provides the necessary basis for the estimation of the curvature 
tolerances. 


§3. THE TOLERANCE IN THE AXIAL THICKNESS OF A COMPONENT 
OR AIR SPACE 


The thickness of a lens is a quantity which is much more difficult to control 
in production than the curvatures of its surfaces, so that the question of tolerances 
becomes very important. For any component the transfer coeflicient, 0A’/dd, 
for its rear surface specifies the change of the aberration per unit change of axial 
thickness. In accordance with the signs we have used, a negative d-change 
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increases the thickness of the component. Formally, then, we introduce an) 


axial thickness coefficient, 0A’/dt, defined by | 
eateag: (sz) es (100) |} 
Ot, dd, | 
the subscript A referring to the general component and the subscript 2 denoting | 
the second, or rear, surface of that component. Thus, without further compu-| 
tation, a set of coefficients is available for the estimation of the tolerances for the 
thicknesses of the components and air spaces. We make use of them by assigning | 
to each component and air space a permissible thickness variation such that the |} 
sum over all the components of the effects on the residual aberrations due to, 
such variations shall not exceed some small fraction of the measures of the residual 
aberrations, 


§4. THE TOLERANCES FOR THE REFRACTIVE INDEX AND 
DISPERSION OF THE GLASSES | 
It was shown in Part I that for any monochromatic aberration, A’, there }} 
may be calculated a differential coefficient, 0A’/0N,, which measures the change’ }} 
of this aberration in the image formed by the system per unit change of refractive } 
index of the component, h, of the system. These coefficients provide full infor-f} 
mation as to the limits within which the refractive index of each component must 
be controlled for satisfactory production of the system, and permit a system off} 
tolerances to be established. As before, we assign a permissible variation of |] 
the refractive index to each component such that the sum over all the components}f} 
of the effects on the residual aberrations due to the permitted variations is less | | 
than a prescribed small fraction of these residuals. As regards the variatiol| 
of the dispersions of the glasses, the coefficients, 0Lch’/0P, and dT ch'/OP,,\I| 
developed in Part I furnish the basis for tolerances, but generally speaking, the] 
variation of the dispersion of a glass between successive melts is negligible, sof 
that these tolerances are seldom required. Frequently a more serious effect of] | 
the variation of refractive index is the resulting change in focal length of the 
system. ‘This is dealt with in the next section. 


§5. THE TOLERANCES FOR THE CONTROL OF THE 
FOCAL LENGTH OF THE SYSTEM 


It is often important to ensure that the focal length of a particular system is}] 
controlled within fine limits during production. For this we must know the] 


effect of the variation of curvature, thickness, and refractive index on the focal | | 
length. From the results of the paraxial ray trace 


f =y|ux 
from a ray incident parallel to the axis at a height y on the first surface. Hence 


of’ 12 Ou,’ f 
dc, (y/Mx Te || 


~ (f' uy’) E. -seeee(101) 


C; 1} 


| 


Transfer coefficients for use in the design of lens systems: IV 429 


Corresponding expressions hold for df’/dn; and @f’/dd,. For the effect of the 
variation in the thickness of a component we have 


o = (f' /uy (Fe) aus (102) 


and the effect of eee variation is given by 
Oh eae NG a Caen! 
ON, Ci NG On Ne 
The transfer coefficients in equations (101) to (103) are quickly computed as the 
values of Ou,’ /0c,, etc., have been calculated in the general computation for the 


paraxial ray. ‘Thus a series of tolerances may be set quite simply for the control 
of the focal length of the system. 


§6. THE USE OF COMPONENTS OUTSIDE THE SPECIFIED TOLERANCES 


When a large number of optical systems of a certain type are being produced, 
a certain percentage of rejected components with dimensions outside the specified 
tolerances is inevitable. The intelligent use of the information provided by the 
transfer coefficients calculated during the design enables some at least of these 
rejects to be sorted into sets of optics which will give a satisfactory performance. 
Components rejected on account of non-spherical figure are excluded from our 
considerations. 

The rejected components are classified according to the type and amount 
of their departure from specified dimensions. For each aberration a limit is 
set for the additional amount of this aberration which can be admitted on account 
of departure from exact specifications. ‘This will be some small fraction of the 
calculated residual aberration. A table is prepared which sets out the aberration 
changes for departure from the specified dimensions in terms of the units used in 
the workroom. ‘The trained worker then derives from this table the amount 
of each aberration introduced into the system if a certain reject component, A, 
is employed. A study of the transfer coefficients in the table will now reveal 
whether some departure from dimensions of another component, B, will introduce 
amounts of the various aberrations which will compensate those introduced by A, 
and thus bring their totals within the limits prescribed. If this is possible, a 
reject of type A is paired with B, etc., and in this. way use may be made of the 
rejected components. This kind of attempt to salvage components is probably 
made in every optical shop by trial-and-error methods or by bench tests. ‘T’he 
use of a table of coefficients as suggested organizes such attempts intelligently, 
giving a reliable guide to what is possible in this respect. 
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A SYSTEM OF TRANSFER COEFFICIENTS FOR 

USE IN THE DESIGN OF LENS SYSTEMS: 

V. TRANSFER COEFFICIENTS FOR THE ASTIG- 

MATISM AT SMALL APERTURE AND FINITE 
OBLIQUITY 


By Fo eD, CRUICKSEDANEG 


University of ‘Tasmania 


MS. received 26 February 1945 


ABSTRACT. ‘Transfer coefficients are developed which specify the changes of the | 


positions of the tangential and sagittal foci in narrow pencils of any obliquity due to small | 


changes in curvature, refractive index, and axial separation which may be made within |} 


the system. 
§1. INTRODUCTPON 


developed for the tangential aberrations only of a lens system. ‘To com- 


te the preceding papers of this series, transfer coefficients have been ]} 


plete the development it is necessary to establish a means of handling the | 


astigmatism in a similar manner. ‘The full solution of this problem, however, 
requires the extension of the foregoing methods to the general skew trace. Not 
having had opportunity to give much attention to this problem as yet, the writer 
has used a simple analysis of the astigmatism at small aperture in pencils of any 


obliquity, checking the astigmatism at full aperture from time to time by the usual |f} 
skew trace. ‘This constitutes quite a useful designing tool, as the changes in the |f) 


astigmatism at small aperture due to changes within the system can be estimated 


for the differential correction process. The transfer coefficients for the shifts 
of the foci of narrow fans of tangential and sagittal rays are deduced, and this 


provides a basis for estimating the astigmatism changes produced by any altera- 
tion made within the system. 


§2. THE TRANSFER COEFFICIENTS FOR THE SHIFT OF THE 
TANGENTIAL FOCUS 


In addition to the principal ray of the oblique pencil of given obliquity, we 


trace another ray of the same obliquity close to the principal ray and calculate ]} 
for it the single surface coefficients and the transfer coefficients in the normal | 
way. We denote this ray by the subscript cp, the letters standing for close | 
principal ray. From the trace of the principal ray we have, as in Part I, ]| 


equation (45 a), 
en GOL 
— : 
“Op, | Op, 


| 
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The two rays pr and cp will emerge from the last surface of the system at different 
points, and it is obvious that we may specify the separation of the emergence 
points by 

OD preemie (Cee ey e—nty OC ammen PEED Tee (104) 

If a change is now made at some surface within the system, the two traced rays 
will emerge from points with a different separation and will intersect in the final 
image-space under a new angle. Thus if we can calculate the changes in 69’,,; 
and 5U;,’ resulting from the alteration made within the system, we can calculate 
the new value of ¢,,’ specifying the position of the tangential focus. 

Let us suppose that a small curvature change is to be made at surface z within the 
system. Differentiating equation (104), 


7) ‘ LO ae ; Ol a : 
ac. (SP orn) a Li, 3c, (Ven = Unie 45 3c, (Ue iad Celie: 


For the term on the left-hand side of this equation we have 


0 , a OP ink en OP prk 
ae rk) a ac; dc; ) 


and combining these two equations we obtain 

Ot,” 1 Pe k Op, k , OU core 0 rat) 

— =—— pe SP) it, | Ke ety]. Cw... 105 

Goma no Un: [( ac, 0c; EA\toc: dc; ey 
The corresponding expressions for thickness changes and refractive-index changes 
may be written down, giving 


ee ep oe ave (? Ucok _ ae) | aly (106) 
Dome SUL \wode > “Od, Gia eed 

Ot; hic 1 Penk Eee) jae (- Uenk = oe) | e(L07 
On OZ i( On; On; *\ an; Gna te ey) 


It will thus be seen that the transfer coefficients for the ¢-focus require little 
calculation beyond the work involved in the trace and coefficients of the ray cp. 
If changes are made during the process of differential correction which are too 
large to be treated as differentials, an improvement in accuracy may be obtained 
by calculating the new values of dp,’ and 5U;,,’ resulting from the change and 
Jeducing the new value of ¢,,’ directly from them. ‘Thus, for a curvature change 


c;, we have 


poison 

ro See (“pe 2 “Boek Sc ves ve (108) 
Zn EON a 

new (8U;/)=8U;,! + (Se a oe] a Wee (109) 


ind the change in the ¢-focus is given by 


1 OP x new / 10 
Of a= (Pi: now pC ls aN el Rae one (110) 
‘ (6 U;, Nee ‘ 


i 
| 
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§3. THE EFFECT OF CHANGES AT A SINGLE SPHERICAL SURFACE) 


ON THE FOCUS OF A NARROW FAN OF SAGITTAL RAYS | 
The position of the sagittal focus behind the successive surfaces of a system 
is easily calculated by the well-known s-trace which is associated with the ordinary 


trace of the principal ray of the oblique pencil under consideration. We propose | 
to use the normal relations of the s-trace to develop the transfer coefficients for | 
the sagittal focus. The position of the sagittal focus behind any surface of the] 


system is given by 


Aa esi cose ae (111) 
s r S 
If we now write c= lis,sor= Lis, 5 og ee (112) 
and G=(N' cost’ =N cost)iN = 5 ere (113) 
equation (111) becomes 
o-=Gotne:; 9) Va eee (114) 


‘This procedure clears the basic equation of reciprocals and makes it easier’ | 


to handle in the subsequent differentiations. Let us suppose now that a small 


curvature change, dc,, is to be made at the surface 7 of the system. The Sly CH 


on the s-focus behind the Suge will be determined Ne the derivative | 


The change o; due to a Sail Paes: fir anes is of second order only, so}f| 
that we may neglect the last term on the right-hand side of the above equation, and} 


hence 


To be able to compute the derivative da;'/dc, we require an expression for 0G;/dc. 
This may be obtained by differentiation of equation (113) and the refraction}! 


equation NsinJ=N’sin I’. Thus, 


0G ol’ 
Afi —sin] ap tusind 


: cos I 
= —sinl’n 
Cc 
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Differentiation of the computing equation 
(L—r)sin U=rsinI 


gives S = Lsin U/cos J, 
and hence we have 
aG _aGal 
dc Ol a 
ee? 0U’\ LsinU 
= sin’ (- 1s) cos I 
rc) Uke 
= sind’ ore em SE he Bo (117) 
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The computation of 0G/dc at any surface thus involves only one operation beyond 
the calculation of the ordinary single surface coefficients, which will have been 
computed already for the principal ray. Two further operations complete the 
calculation of do’ /dc. 

To determine the effect of a small change in refractive index at any surface i 
of the system, differentiation of equation (114) gives 


do; 0G 
Be ee svelecaie sé (118) 
in which, by differentiation of (113), we have 
2 = —sin I’ — —cosI 
On 
=-—siny’ oh —cos I 
os I 
= Spa. eee Shee (119) 


In this case only four operations are required to calculate do’/dn for any surface 
after the single surface coefficients. 

Finally, if a small change is made in the axial separation of two consecutive 
surfaces of the system, differentiation of equation (114) gives 


do; 0G: do, 

ad, = Ca aus eee (120) 
Expressions for the two derivatives on the right-hand side of this equation may 
be deduced quite easily. Thus, differentiating equation (113), 


Uc OL: ol 
ane ~sin =; + usind 
rane t gieeols 
= sin ad an 
al’ Or noose ol 
eae ad cos I’ ad 
aa et dU'\ al op 
= sinl (1 abe 
= sin/’ sin pe takes (121) 


op 
To obtain an expression for the second derivative we differentiate the transfer 
formula of the s-trace, which is 


Pet D tT) ge Oth a Ge on (122) 
do; 1 Os; 
and thus ad. es 2 a 
oD; 
— 2 Mi 
ria eheod: 
Bev see u;| 1 eee. send (123) 
: : cos I 


The last step follows simply from the well known relation between D,_, and d;. 
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Often the bracketed expression may be replaced by unity with a good degree ofl] 
approximation, and is accurately unity for a plane surface. Equations (120), ! 
(121) and (123) permit the calculation of do’/dd at any surface of the system. | 


§4. THE TRANSFER COEFFICIENTS FOR THE SAGITTAL FOCUS 
OF A NARROW PENCIL 


It now remains to develop a method of computing the transfer coefficient 
for o’-changes so that the effect of a change within the system on the final value of} 
a,’ may be determined. Since G is a function of the angles of incidence and] 
refraction, its value at any surface will change with change of the point off] 
incidence. It will be found, however, that in general the change of G along the] 
path of the pencil due to small alterations within the system has a very small} 
effect on the transfer coefficients, and, on account of the simplification obtained} 
by neglecting this factor, it is proposed to develop a set of transfer coefficients ]} 
which leave it out of consideration altogether. 1 

Suppose that by some means a change do,’ is produced in the narrow pencil} 
refracted at surface 7. To determine the effect of this change at the following} 
surface of the system we require an expression for the derivative 00;+,/00;., | 
We may write the transfer formula in (122) in the form 


0144 = 0; (1 — Dio; 5 ee ee (124) 
and hence ou rt = (1 Dia, Pa ie ea eae (125) 
oe, 


The change which will result behind surface (¢+1) must now be determined. 
From equation (114) we have 


0a; 
ae Se 
ae Ness a ee ee (126) 
Oy41 
da; 0a; 
trl i+1 
whence = = Nita; 
do; do, 


and extending this expression surface by surface we obtain 
Oo;  Oo;44 00; 49 0o;, 
36 ar tt ate Seep Fee, 
and the relation between these transfer coefficients at successive surfaces is, 
therefore, 
Oo% _ O41 Oj 128 
dof de, hide a, (128) | 
We now have the complete framework for the computation of the transfer coeffi- |) 
cients. Equation (125) is used first to calculate dc,,,/00,’ for each surface except |f) 
the last. ‘Then at the last surface 00,’ /dc,,’=1, and equation (128) is used to 
calculate 0o,,'/0c;,'_,, and the computation may be continued surface by surface ]} 
through the system. Finally these basic transfer coefficients are used in conjunc- |) 


tion with the single surface coefficients, giving expressions of the form 


Oo, _ Aoz,’ Ce : 


dc; 00, Oc, 

and returning to the quantity s’ instead of its reciprocal, we have 
OS; 19, Oc, 0a,’ ‘ 
BGs * Gof ioaeeea See (129) 


with corresponding expressions for 0s, /dd; and Qs,’ /On,. 
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§5. A NUMERICAL EXAMPLE 
In the course of correcting a certain photographic objective having ten 
surfaces, a change of curvature of 0-001 mm.-! was made at the seventh surface, 
the radius of curvature of that surface changing from 1091 mm. to 521-78 mm. 
In a pencil of approximately 17° obliquity the following results were obtained: 


OSi. Por Wis, LE 5A 
OC, OC, 
and hence 
ds, = — 2-113, at,’ = — 2-416, 
giving 
news; = 47:273, newt,’ = 47-661. 
Using the alternative relations (108) to (110) there was obtained 
diy, = 2°33, 


newt,’ = 47-742. 
A full trace of the altered system gave the final positions of the two foci as 
Sy = 947-338; t, = 47-801, 
which represents an accuracy of about 3°%. 


A TRANSFER METHOD FOR DERIVING THE 

EFFECT ON THE IMAGE FORMED: BY AN 

OPPICAL SYSTEM- FROM RAY | CHANGES 
PRODUCED AT “A GIVEN SURFACE 
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ABSTRACT. A transfer method is described depending on the invariance of a function 
of the position change of successive foci of a narrow pencil as it passes through an optical 
system. 

The method described is approximate, as the invariant is applied to wide pencils, but 
by reason of its speed it is a useful tool in design. 


§1. INTRODUCTION 

N a paper by M‘Aulay and Cruickshank (1945) a method is described by 
| which the effect on the final image formed by an optical system may be 

computed when alterations are made at any surface of the system. The 
method naturally breaks into two parts. In the first part the changes which 
rays suffer as they pass through the surface under consideration are evaluated, 
and in the second part these changes are transferred to the final image formed 
by the system. 
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The present paper describes another method of transferring the changes} 
produced at any surface to the final image. It is based on a theorem of a form| 
similar to Lagrange’s, but is applicable to pencils of any obliquity.* If the} 
theorem is applied to narrow pencils it is exact for differential changes, but inj} 
the use here described it is applied to wide pencils and is approximate. It has} 
been used principally as a primary computing method where speed is more} 
important than exactness, and it has been found of considerable practical value.}} 
It is believed that its potentialities are greater than the particular development | 


discussed. 


in position due to an alteration made in the constants of one surface of the syster 
through which it is passing. Figure 1 shows a narrow pencil ABC focusing} 


BY’ Z 1s = y 
Jae ai 


Figure 1. 


the ith surface causes the pencil to focus instead at B.* The co-ordinates of thd 
change are w, and z; measured along and at right angles to the original axis of the 
pencil. They are given the positive sign when they have the direction shown} 
in the figure, as this relates them conveniently to the usual sign of U, the angle} 
under an incident ray. If y; is the angle between the extreme rays of the pencil} 
and WN; is the refractive index of the medium between the 7th and the jth surfaces | 
it may be shown that | 


(1) Nyz remains unchanged as the pencil is followed through successivd} 
surfaces in its original and changed form, and 


(ii) Ny?w also remains unchanged. 


The second relation holds, however, only if the change is such that a spherical, 
wave remains spherical after the change. It is a useful relation for some purposes i 
but will not be further considered here. 

* 'This theorem, given in § 2, was discovered by the author. 'The referee has pointed out the | 


it was known independently (see, e.g., v. Rohr, Formation of Ima | 7 
; evs < ges in Optical Instr 
1920, p. 174), and no proof is therefore given here, p om 
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It is important to note that the change of a ray is not being followed. The 
%8 are co-ordinates of the changes in the successive focal points of a pencil, but 
the actual rays which have been followed in the above argument change from 
surface to surface. Unless this point is kept clearly in mind it is easy to make 
serious mistakes when the principle is used as a computing tool. 


Alternative statement of the principle. Figure 2 shows a pencil with its axis 


.Figure 2. 


displaced a small distance 2’. It is evident that if «’ is the angle through which 
the axis is displaced, 


N’y' , — Ne = N'a«'dp", 


where S’ is the length of the axis of the pencil from surface to focus, 5p’ is the 
width of the pencil at the surface and «’ the angle turned through by the axis 
of the pencil due to the alteration. ‘This form is often best suited for obtaining 
he value of the invariant at an altered surface. 


SoA PETC ISLON OBS THE APRINCIE TE EMO. EE, aRANS FER Oe 
DIFFERENTIAL CHANGES PRODUCED AT THE iTH SURFACE 
OF A SYSTEM 


In the paper by M‘Aulay and Cruickshank quoted above, a change in the 
iberrations of an optical system due to an alteration of the 7th surface was dis- 
ussed. ‘The system is primarily analysed by considering six rays, three 
ixial (marginal, zonal and paraxial) and three oblique (a, pr and 5). ‘The same 
ays are considered in the method to be described, but they are treated as the 
xes of six pencils. It is shown in the paper quoted that the change in angle 
inder a ray refracted at a spherical surface can be given as a function of the 
urvature or refractive-index change. In what follows, it will be assumed that 
uch a surface alteration has produced changes dU’ in the angles under the refracted 
ays, and that these have been computed. 
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The transfer theorem is used to convert these dU to displacements 2’ per- 
pendicular to the directions of the changed rays in the neighbourhood of the final | 
image plane. It is evident that 2’ is very nearly the co-ordinate of the displace- }} 
ment of the focus of the narrow pencil which in the notation of the last section |} 
would have been written z;,’, where k is the final surface. : 


Approximate evaluation of 2’. The theorem is used in the form 
N’a’dp’ is invariant through the system. i] 


This invariant is found for the 7th surface and used to obtain 2 after refraction’ | 
at the last surface by the relation 


/ / VA / / / red y 
N,' 4; 8p,’ = Ny! %' dp,’ = Ni, = Pre 
3p; 


or z = N, a8 5! PyHe 

where the refractive index after the last surface considered is unity. 
N,' is available from the original data, «’ is the change dU; which has been ]} 

supposed computed, and S;,’, the distance measured along the ray from the last }} 

surface to the final image, can easily be obtained from the trace. It remains to) 

evaluate 5p,'/8p;. . 1 
It is here that the most serious approximation is made. The object is to find 


sufficiently good values of 8p’ from the basic six-ray trace without much extralf] 


Figure 3. 


computation. In the tracing method used, the value of L’ sin U’ =p’ is, inci- | 
dentally, found for each ray. his is the perpendicular distance from the pole. 
of a surface to the ray under consideration. 

Figure 3 (a) shows marginal, zonal and paraxial rays; figure 3(b), a, pr and 
b rays and the perpendiculars dropped from the pole on them. The following} 
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values are taken for (5p'/3p,’): 
Px Pick = Px |Prn 
Op. [SP:x = Pur |Purm 
OP Opin = 2( Pu’ —Pa’)|(Par — Pa )n— Pa [Pow 
Pprl Pore = (Pa —Po')/(Pa' — Po’) 
Sa [Par = (Pa — Pr) (Pa — Por) 
+ 3[(Pa’ = Ppr)l(Pa' — Por) — (Por — Po’ )I(Por—Po als 
OP» Pox = (Por —Po')I( Por —Po dx 
— 31(Pa' — Ppr)|(Pa' —Por)e—(Ppr Po’) Por Po )al- 
{t is evident that other extrapolations for the marginal, a and b rays might be 
used. ‘I’he rather crude forms given are easy to compute and reasonably good. 


§4. THE ABERRATION CHANGES FROM THE 2 


Spherical aberration and curvature. ‘The aberration changes can be stated 
simply in terms of 2’. ‘The changes in spherical aberration and tangential 
curvature of field will be written down without further discussion as their 
derivation follows the same lines as that given by Cruickshank, Part I (1945). 

For spherical aberration 

aly =dl —dLy =2, [u —&y coséec Uy, , 
with an exactly similar expression for zonal spherical aberration. 
For curvature 
dX 7 =(2, cos U,’ —2,' cos U,') cosec(U,' — U,') —2,/ cosec U,’. 

The notation is that of Conrady’s Applied Optics and Optical Design or the 
papers quoted above. A change of sign, however, has been introduced for X7’, _ 
and the signs for coma and distortion will also be changed. It has been found 
convenient in this work to adopt a positive sign for an aberration that arises when 
parallel pencils of small inclination to the principal axis pass from air through a 
convex glass surface. 

Coma and distortion. ‘The co-ordinate changes perpendicular to the principal 
rays can be used direct to estimate these aberrations. ‘The change in coma 1s 
sufficiently well given by the change in distance between the principal ray and 
the median of the a and 6 rays. It is given by 

d(coma’) = 2,,,—3(2a— 2, )sec 3(U, — U,’). 

The treatment of change of distortion, involving as it does an imaginary ray 
not included in the trace, requires a little more discussion. An ideal ray is 
considered which is identical with the principal ray till it strikes the surface at 
which an alteration takes place. Subsequently this ray is treated as though it 
were refracted according to the paraxial form of the computing equations. The 
change in angle under the ideal ray after refraction at the altered surface therefore 
differs for that obtained for the principal ray. This change of angle is 
computed and is transferred to a 2’ change in the neighbourhood of the final 
image, which is written 2),,, and is given by 


Sone =r IN Kno (Pa Pon) ee 
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The change in distortion is now given in the same manner as the change in 
coma by 


| 
| 
| 
| 
: 
| 


A( dist) =2 55-26 
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ABSTRACT. Values are calculated for the brightness of a total radiator at temperatures 
between 2000° Kk. and 3120° k., using the standard C.I.E. relative luminosity curve and 
the Planck radiation law. A simple logarithmic equation is found to express the relation 
between brightness and temperature to within +1% over this temperature range. The 
ratio of the new international candle to the lightwatt is calculated from the results. 


§1. INTRODUCTION 


radiator as a function of temperature. From the literature it appears 

that nothing has been published on the subject since the adoption of the 
C.L.E. relative luminosity curve in 1924. Hyde, Forsythe and Cady (1919) pub- 
lished some values of the brightness for the range 1700° k. to 2650°K.; they used 
their own values for the relative luminosity function, however, and these differ 
from the standard values later adopted. ‘hey expressed the results in terms of 
a rather complex four-term equation. Nernst (1906) expressed the relation 
experimentally found between the brightness (B) and temperature (T) of a total 
radiator for the range 1460° k. to 2280° k. by the following equation :— 

log B= —A/T+C, 

but the relative luminosity curve for the observer concerned is not recorded. 

Values for the brightness of a total radiator have been calculated from the 
C.I.E. relative luminosity data and the Planck equation for fifteen temperatures 
ranging from 2000° x. to3120°K. The following approximate formula, which 


has the same form as Nernst’s relation, has been obtained for the brightness of a 
total radiator over this temperature-range :— 


lop), B= 4:275 — 10860) (1) 


where B is the brightness of the total radiator (lightwatts em.-*) and 7's its 
temperature (°K.). 


[ is sometimes desirable to use an expression for the brightness of a total 


Relation between brightness and temperature of a total radiator AAI 


§2. CALCULATION OF BRIGHTNESS 
The brightness of a total radiator at a given temperature is calculated from the 
equation 
1 Pic.) 
B=a- | Vd aoe ee Re (2) 
0 


7 
where B is the brightness (lightwatts cm.’). J,dd is the total energy flux 
per unit area from a total radiator (watts cm. *) between the wave-lengths A 1o 
(microns), and V, is the value of the relative luminosity factor at the wave- 
length A, as given by the 1924 C.I.E. standard luminosity curve (Judd, 1931), 


1 é ; : 
The factor — appears because the brightness of an element of surface, i.e. 
aT 
the ratio of intensity to projected area in a given direction, bears the ratio — 
7 


to the total luminous flux per unit area from the same element of surface (assuming 


that the surface has the properties of a perfect diffuser). In the present case the 
Pie.) 


expression | J,V,dX represents the total luminous flux per unit area of surface, 
0 
since it is derived from J,dA, which represents the total energy flux per unit area 
; oe 1S 
for a small wave-length range. Hence the brightness is given by — | J,V,da. 
MA 


The values of J, dA are derived from Planck’s formula, 


C5 dA 


J,dx = got _ | Se Be ee iw erersilers (3) 


where C,, and C, are constants and J, dA, X, and T have been previously defined. 
Convenient tables for deducing values of B from equation (2) are those 
of Skogland (1929). These tables employ the values for C, and C, given in 
the International Critical Tables, namely C,=3-703 x10 erg sec>1 cm*?, 
C,= 14330 micron degrees. ‘The range of temperature considered is 2000° kK. 
to 3120° x. ; this is covered in steps of 20° in table 1 and 80° in table 2. Table 1 
gives at each temperature, besides values of J,/Jo.5) for a series of wave-lengths, 
values of Jo-59/C, to 5 significant figures (units: microns °). From these latter 
values, taking C, =3-703 x 10> erg sec>1 cm??? and V».5,=0-7570 (Judd, 1931), 
we calculate the value of Jo.59Vo.59 at a series of temperatures, at intervals of 80°, 
from 2000° x. to 3120°K. Skogland’s table 2 gives at each of these temperatures, 
for a series of wave-lengths from 0-401 to 0-76 in steps of 0-01 p, values for 
J,V,|(JV) max, Where (JV)max is the maximum value of J,V, at the given tem- 
perature. Reading off the values for Jo.59Vo59/(JV max, and substituting the 
value already obtained for Jo.59 Vo.59, we find (J Vmax for the given temperature. 
Table 2 gives also the sum of the values of J,V;/(JV )max, at each temperature, 1.e., 
1=0°76 


~y aly V/(J eae, 


A=0°40 uw 


the summation referring to steps of 0-01 u. ‘The values at the wave-length limits 


442 E. F. Caldin 


are not more than about 0-03°% of the maximum values, and the sum does not 

differ by more than 1 in 104 from 
A= 

» d,, V,/(J UA ee | 

2=0 | 

Multiplying the sum by 0-01, we obtain, therefore, a value (in microns) of} 

i: . [oie | 

[J,V,\(IV max] @, ive. arama TV 

I, A 4l( ma ] Chae . Ameee 

which is accurate at least within 1 in 104. Multiplying this by the value of 


(JV) max already found, we obtain the value of | J,V,dXd; and hence, dividing 
JG | 


by =, the brightness of a total radiator, at the given temperature. 


§3. RESULTS AND DISCUSSION | 

In the table below are given values for the brightness B, in lightwatts per 
sq. cm., and for log,,) B, at a series of temperatures from 2000" k. to 3120°K. 
The value of log,, 2 is found to vary nearly linearly with 1/7, and the following 
formula has been tested :— | 
logy) B= 42750 —10360-0/T, 7) eee (1) 

Values derived from this formula are given in the table and compared with] 
the values directly calculated. | 


log 1B 

Temperature| Brightness B log,)B calculated Difference 

(Giza) (lightwatts cm>?) (a) from egn. (1) (a)—(b) 
(b) 

2000 0:0706 2-8490 2:8450 + 06-0040 
2080 0:1137 1-0556 1:0538 +0-0018 
2160 0-1767 1:2473 1:2472 +0-0001 
2240 0-2665 1:4257 1-4268 —0-0011 
2320 0-3910 1:5922 1:5940 —0-0018 
2400 0:5595 1:7478 71-7500 —(-0022 
2480 0-7827 T-8936 1-8959 —0-0023 
2560 1-0730 0-0306 0-0328 —0-0022 
2640 1-443 0-1594 0-1613 —0-0019 
2720 1-910 0-2809 0:2823 —0-0014 
2800 2-486 0-3956 0-3964 —0-0008 
2880 Selgy 0-5040 0-5042 —0-0002 
2960 4-040 0-6068 0-6061 +-0:-0007 
3040 5-060 0-7042 0-7026 -+-0-0016 
3120 6°262 0-7967 0-7942 +0-0025 


__ Itis evident that the values of B are given within + 1% (an accuracy sufficient, 
for many purposes) by the empirical equation (1). 

The effect on the formula (1) of using other values for C, and C, may be} 
noted. A change in C, produces a proportionate change in all the values of 


J di, and hence in | J, Vdd, and so can be met by changing the constant 4:275 to 
0 
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[4-275 + logy) (C,/3-703)], where C, is expressed in erg sec. em: 4~< 10>; 
A change in C, affects the temperature-dependent term; since exp (C,/AT)>1 
for the ranges of temperature and wave-length concerned, it is only necessary 
to change the constant 10860 to (10860 x C,/14330), where C, is expressed in 
micron degrees. Thus if we adopt the recently suggested values of Birge (1941), 
namely C,=3-7430 x 10-> ergsec.-'cm.t? and C,=14384-8 micron degrees, 


the relation becomes 
log, B=4-279 — 10902/T. 


$4. PHOTOMETRIC UNITS: THE RATIO OF THE NEW 
INTERNATIONAL CANDLE TO THE LIGHTWATT 

It is possible to convert the above values of brightness directly from lightwatts 
per sq. cm. to “new international candles”’ per sq. cm., and to calculate the ratio 
of the new candle to the lightwatt, because the new candle is defined in terms of 
the intensity of a total radiator at a fixed temperature. This is of interest because 
the new candle is.likely to come into force in the near future as the international 
unit of luminous intensity. By definition, the luminous intensity of a total 
radiator at the temperature of solidification of platinum is 60 new candles (Inter- 
national Lighting Vocabulary, C.I.E., 1938). The brightness in lightwatts 
per cm? of a total radiator at this temperature (taken here as 2045-9° x. 
(Wensel, 1937) is found, by interpolation from the directly calculated figures in the 
table above, to be 0-:0932. ‘The ratio of the new candle to the lightwatt is thus 
60/0-0932 = 643 for the values of C, and C, used in computing the figures in the 
table ; and to convert the values given above for the brightness of a total radiator 
to new candles per cm? it is only necessary to multiply each value by 643, or to 
add 2-808 to its logarithm. 
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Note added in proof. Mention should be made of the work of Ives (1926), 
who calculated values of the brightness of a total radiator at temperatures from 
1200° x. to 10000° x. In these calculations, an approximate formula was used 
0 express the C.I.E. relative-lurninosity data, and errors of the order of 
[°%, are likely. The author was unaware of this work when the above paper 
vas written. He has also heen informed of a paper by W. Geiss (Licht, 
[2, 33 (1943)), but has not been able to see a copy. 
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REVIEWS OF BOOKS 


Waveform Analysis, by R. G. Mantey. Pp. viit+275 and 3 plates. (London: 
Chapman and Hall, Ltd., 1945.) 21s." 


his book is not a treatise on harmonic analysis in the conventional sense but, as the | 
author correctly emphasizes in the sub-title and preface, a practical guide to the rapid | 
interpretation of wave records. Its primary aim is to explain, against a background of | 


more formal analysis, methods, involving merely inspection and simple measurement, 


for estimating the frequency, amplitude and phase of the principal harmonic components _ 
of a large class of complex wave-forms which occur especially in vibration engineering, | 
‘There are ten chapters, entitled : I, Sine waves in combination ; I1, General properties of | 
harmonic series ; 11, Basic analysis of recorded waveforms ; IV, The envelope method ; | 


V, Method of superposition ; V1, Fourier series: mathematical analysis ; VII, Numerical 


methods ; VIII, Mechanical and other aids to analysis ; 1X, Practical requirements for } 


waveforms ; X, Lissajous figures. 
Chapter I illustrates graphically and trigonometrically the synthesis of wave patterns 


(especially beats) from two or more sine components, and, together with Chapter II, 


provides that foundation of experience of the anatomy of wave-forms which is necessary 
for a ready diagnosis of harmonic content. 

Chapter III seems a little dilute, and might be made more concise with advantage, 
but it serves to pave the way to the actual measuring of records and hence to Chapter IV, 
which forms the pivot of the book. This chapter, essentially diagnostic in character, 
shows how the construction of ‘‘ envelopes ’’, to touch the crests and troughs (or, more 


generally, to pass through any repeated characteristic kinks) of the undulations, makes it |] 


possible to identify the harmonics and to determine their characteristics. A lucid explana- 
tion is given of the surprising amount of information which may be extracted from a 
wave-form by an experienced interpreter without recourse to elaborate analysis. The 
method appears to have especial application to the routine examination of vibrograph, 
strain-gauge and similar records which are nowadays obtained in certain industries in 
numbers too great to permit of detailed analysis. A systematic procedure for applying 
the envelope method is developed in a number of tables somewhat reminiscent of chemical 
analysis. 

For certain frequency ratios amongst the harmonics (such as 3:2: 1), the envelope 
method fails. Chapter V shows how it is possible to meet these cases by dividing the 
complete cycle of the given wave into a number of equal parts which, by suitable summation 
and differencing, permit partial, and in favourable circumstances complete, separation of 
the harmonics. 

It may be remarked that the envelope technique and similar simple methods, while of 
great service in resolving a wave composed of relatively few harmonics of comparable 
amplitude, do not promise so much help to the physicist, meteorologist and economist, 


who will still have to resort to Fourier, periodogram and correlogram analysis when con- | 


fronted, as they frequently are, by waves having a long retinue of progressively diminishing 
harmonics. 

The theoretical Chapter VI treats the formal representation of functions (including 
the usual examples of square and saw-tooth profiles) by Fourier series, and gives an adequate 
sketch of the mathematical background. While very properly leaving full rigour to the 
mathematicians, it does somewhat more than merely give plausible derivations of formulae. 
One cannot help regretting, however, that the nature of the book scarcely permits pre- 


sentation of the connection between the author’s ‘‘ periodic existence functions ”? and their _ 


analogues, the Cauchy discontinuous functions. Gibbs’ phenomenon at discontinuities 
is treated briefly, and Fourier integrals are mentioned but not applied. 


Chapter VII gives a clear account of the numerical Fourier analysis of waves specified | 
by equidistant ordinates, including computation schedules for schemes using 24 and | 


48 ordinates. 
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Chapter VIII describes in detail a versatile drawing-table type of mechanical harmonic 
analyser, and briefly mentions alternative devices, including electrical filter circuits. No 
attempt is made to survey the very rich field of mechanical and electrical aids which now 
exists. 

General comments on the practical problems of recording and allowing for response 
distortion are given in Chapter IX, and the book closes with a chapter on the characteristics 
of two vibrations compared by means of Lissajous figures. - 

There are several appendices dealing with mathematical points, a list of elementary 
trigonometrical formulae, a useful table of sines and cosines adapted to facilitate harmonic 
analysis, a glossary of terms, a bibliography and an index. 

A minor criticism of the presentation is that elementary details of mathematical manipu- 
lation are occasionally given more fully than is necessary if the reader is sufficiently mature 
to appreciate the chapter on Fourier analysis ; but this may encourage further reading, 
and it would therefore be beneficial if the bibliography were somewhat enlarged. Viewed 
as a whole, the book is well written and competently composed, and should be of con- 
siderable value to those concerned with vibration engineering. There are few misprints, 
of which the most obtrusive is the insertion of the terminal “‘s”’ of Lissajous in L’ Hopital. 
The printing and binding are quite good. M. S. JONES, 


Introduction a étude des champs physiques, by J. GRANIER. Pp. iv+251. 
(Paris: Dunod, 1941.) 


This book is not an introduction to field theories as now understood, but a general 
treatment of steady-state boundary problems where the differential equation is that of 
Laplace or of Poisson, that is of problems in electrostatics, magnetism, heat conduction, 
capillarity or hydrodynamics. ‘The treatment is relatively elementary, and is confined 
in the main to problems in two dimensions, or those with spherical symmetry, so that 
there is no call to introduce Legendre functions, spherical harmonics or Bessel functions. 
‘The French nation has given the world some of its greatest geometers, and the tradition 
of geometrical reasoning has always been strong there, so that it is not surprising to note 
that, in this book, transformations from one co-ordinate system to another are seen from a 
rather more directly geometric aspect than they usually are in books written by Englishmen, 
who tend to see them as algebraic transformations. 

After a preliminary chapter in which conjugate complex functions and systems of 
orthogonal curves are considered, the book is divided into chapters, each dealing with 
some phenomenon in which Laplacian fields present themselves, beginning with electro- 
statics. By taking the potential due to a point charge, and that due to an infinite straight 
line (a point charge in two dimensions), other cases are built up, the most elegant being 
that of a circular cylinder surrounding two line charges of opposite sign, a case which is 
found when considering the case of four line charges, by noting that there is one equi- 
potential which is circular. A few cases are solved by means of conjugate functions, 
but these are generally introduced experimentally “to find out what cases they will 
represent’, rather than deliberately to solve particular problems posed beforehand. 
Even after the theorem of Schwartz on the behaviour of functions at branch points has 
been introduced, little use is made of it, though in fact it enables us to select with certainty 
the transformation appropriate to any plane problem with rectilinear boundaries. ‘The 
chapter on magnetism finds its interest of course in doublets and other multiple charges, 
and does not carry the matter very far. The field due to a sphere in a field which would 
otherwise be uniform is worked out, but the case, so dear to most text-book writers, of an 
ellipsoid in a uniform field, is omitted. The misprints in the middle of page 121, where 
R? is omitted from the denominator of two successive equations, might mislead beginners. 
In the section on hydrodynamics we have such problems as the flow past a cylinder, and 
attention is then turned to cases where there is a free stream-line, as when a jet issues 
through a slit. No use is made here of the method devised by Kirchhoff of introducing 
a subsidiary plane in which the co-ordinates are the absolute magnitude and the direction 
of the velocity, but a great deal is extracted from the simple transformation, attributed to 
Joukowsky, x=3(r-+r-}) cos 8 and y=3(r—r~?) sin 8, which analysts would think of as 
w=(z+1/z). This, with its inverse transformation, and the addition of different fields 
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of flow, yields an astonishing amount of information about flow round obstacles of practical 
interest, such as aerofoils. | 

The last quarter of the book is of a rather different nature, and considers the use of | 
experimental methods for solving complicated problems—measurements of electrical 
resistance of diagrams cut out of metal foil, or tracing potential lines by means of probe | 
measurements in an electrolyte, or plotting out a field by means of iron filings above a | 
magnet, for example. The author has clearly had a very considerable experience with 
these methods, and is able to give advice on experimental arrangements which must be: 
of immense value to any who propose to use the methods. Among them he includes 
the graphical method in which a network of orthogonal curves is drawn freehand, so spaced 
as to divide the field approximately into squares. 

There is no doubt that the book will be of value to many, in showing the essential | 
unity of the diverse branches of physics with which it deals, inenabling them to devise methods. __ 
for attacking problems with which they are confronted, and as a book of reference, giving | 
the solutions of many problems. When we consider that it was published in France in 
1941, we can only imagine the difficulties which must have confronted author and pub- | 
lisher, and must admire them for producing the book at all ; that they produced such a 


good one must almost be a matter of wonder. ji HeAS 


The Quantum Theory of Radiation, by W. He1rLer. Second Edition. Pp. viii+ 
272 (London and Oxford : Sir Humphrey Milford, 1944.) 20s. net. 


Ten years have passed since the first edition of this book appeared, and it has clearly 
established itself in the interval as an authoritative treatment of its extremely difficult 
subject. Although progress, particularly on the experimental side, has been steady, 
there has not yet been any spectacular advance which would call for a revision and 
re-writing of the theory, and in fact the present volume does not differ essentially from 
the first edition, save for interpolations here and there and an addition to the appendix. 

One day, of course, a complete change of approach will be found, which will at one 
stroke give us a quantum theory completely in harmony with the postulates of ‘relativity, 
and which will avoid the difficulties we find at present in connexion with the infinite 
self-energy of the point electron. ‘Till this new formulation is made, it is difficult to 
imagine any more useful, or more instructive, text-book than this, though we may at 


times wish that someone would write a simple text fo serve as an introduction to the 
7. eae 


subject. 
Tables of Elementary Functions, by F. EMpr. Pp. xii+181. (Leipzig: Teubner, 
1940.; photographic reprint published under licence by J. W. Edwards, 


Ann-Arbor, Mich., U.S.A., 1945.) $3.20. 


Sechsstellige trigonometrische Tafeln, by H. BRANDENBURG. Pp. xxiv+304. 
(Leipzig : Lorentz, 1932 ; photographic reprint published under licence 
by J. W. Edwards, Ann Arbor, Mich., U.S.A., 1945.) $5.00. 


It is regretted that, in the reviews of these two books in the July issue of the Proceedings 
(this volume, pp. 368-370), no mention was made of the fact that the sole agent for these 
tables in Great Britain is The Scientific Computing Service, Ltd., of 23 Bedford Square, | 


London W.C. 1. 
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